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Introduction 

Analytical theories such as Brouwer, 1 Vinti, 2 and Kozai, 3 were used extensively for orbit 
predictions in the Space Age early years. The advance in large frame extremely fast computers 
allowed less premium to be placed on computational efficiency, thus special perturbations 
methods e g Cowell etc., could be employed to address the increasingly more stringent accu- 
racy requirements for orbit prediction. However, there has developed an increasing need for au- 
tonomous satellite ephemeris generation by various users, which has been facilitated by the intro- 
duction of mini and micro personal computers (PCs) into this field. This has created a desire tor 
methods which retain computational efficiency while meeting the requirements tor greater accu- 
racy. For example, research developments by Rom, 4 Deprit, and averaging methods by Mc- 
Clain, et al. 6 provide for development of more accurate analytic methods and semi-analytic 

methods which address these needs. . , 

There has always existed a wide area of applications for general perturbation models of mod- 
erate accuracy, providing even more efficient computational speed and storage requirements 
suitable for micro-computers. These needs have been expanded by the increased need to provide 
autonomous satellite orbit propagations on small and slow computer processors onboard a 
spacecraft (with a premium placed on power conservation) or at “on-site” ground locations. 
Thus, there is a continuing demand for moderately accurate general -perturbation models designe 

for speed and low memory size. . 

Hoots’ 7 latest contribution to this genre has noted the legacy of developments in this area. 
Hoots recounts that: “One of the first theories of this type was developed by Hilton and Kuhlman 
in 1966.” The resultant simplified general-perturbation theory obtained from simplification ot 
the work of Kozai for gravitational effects with the drag effect on mean motion taken as linear in 
time is called SGP. “SGP is used throughout the world at radar installations, data collection sta- 
tions, universities, etc., where fast, moderate accuracy satellite predictions are required.” A later 
simplifying theory, called SGP4, was developed in 1970 by Cranford (see Lane and Hoots ), by 
use of Brouwer theory for its gravitational model and a power density function for its atmos- 
pheric model. “SGP4 is currently used by North American Aerospace Defense Command 
(NORAD) for updating and maintenance of the entire inventory of near-Earth satellites. 

Hoots presented a simplified general perturbations model called SGP8 in 1980 which sig- 
nificantly improved prediction accuracy near decay where SGP and SGP4 predictions degrade 
rapidly. The SGP8 theory employs the same atmospheric model as SGP4, incorporating a 
hybrid simplification of Brouwer and Vinti’s general perturbation theories for its gravitational 
model which removes any singularities at the critical inclination. The SGP8 retains most of the 
prediction accuracy except for satellites with larger eccentricity, however transformation from 
osculating to mean elements is not available in closed form (i.e., the theory is not self starting). 
Therefore the seven constants of the theory comprising the drag parameter and mean orbital 
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reference ^ etermined through a leas ‘-square fit of the theory to the first day of a given 

In order to facilitate the efficient use of the general perturbation model presented here on a 
~o P a Ut S; h r P T pre f Se " tS / n ,m P roved algorithm for osculating to mean element conver- 
mafio^methods 2 " matnzant 18 dCVel ° ped t0 emp, ° y With the theory in state estj - 

th - ° f thC sim P llfled g eneral perturbation theories presented including the present offering in 
paper occupy core size of around 8K bytes. In order to demonstrate its adaptation to micro- 
computers the method presented here has been implemented on a very inexpensive PC the 
Timex/Sinclair 2068 color home computer. Standalone BASIC program packages have been 

sxrr s a x ,r elemem d “ — - ** ~ *+£ 

and ^"gt^atio^ propagation UmeriCaI meth ° d C ° nSiderS 34x4 gravitational field 
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Economizing Procedure 

The simplified genera! perturbation theory presented in this paper (hereafter designated "Be”) employs 
implifications of the Brouwer and Lyddane theory for its gravitational model. The theory is further modified to 
nco^orate a computationally efficient algorithm to simulate the drag effects by a retarded linear rate parameter 

o^e The m t n ih mi ' maJ °; aXIS a The P 3 "™!™ derived * observing the semi-major axis decay in mean 
P • . • f eor y _ neglects drag variations within a period, rectifying the constants of the theory at orbit 

maly JmZc. ^ rffca °" ' he mea " a " d ^adra.ic variation in H ,„<? 

Even without the drag algorithm, the method of Brouwer and its revision by Lyddane has been adapted to 
of thlT eCon . on ^ Cal an 5' ytlcal orb,t theory for satellite motion about an oblate planet including J,, L and parts 
the t 4 Z ° na 6 eC f S i e true ar g um ent of latitude was adopted as the fast variable in the theory The choice of 
the tnie argument of latitude as the fast variable simplifies the computation of the oscuIatiM^lin«£nS J 

effects „„ 1 radlaVand c °rT Lyddane " ro " sidered this ada P“->n in relaL to Sr prTma 
tionalconsideradon. a " d '™"““= d ™ accordance with economical compttta 

Lyddane remarks that i" and g" (“mean” mean anomaly and mean argument of perigee) must be used for 

al ^thiTresults in J rdatfvel^h > rg and h^I ^ magn ‘| ude) in his version; however, as demonstrated by Gordon, et 
Thk !! results in a relatively large radial error with respect to Brouwer for moderate values of the eccentricitv 

r ' r: ,h f, 8 ' ,se ) uiar t + ,on ^ i"c^, e 4rSeS 

factors it, the long period variation in the ectSmticity due toTtaJ tte "’“’'i?*! 

Applications 

field^apfdfeations^aMhe fool' of^n anTennVl^ geT" “ °J i,) This P™e S -seful for 

Economized Equations for “Bg” 

nf rh^K-Tf ‘t 16 S,mplifyin 8 Procedures discussed we arrive at the following formulas for efficient computation 
the orbital elements using the same notation and definitions employed by Brouwer. 1 P 


Abbreviated terms: 
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Secular terms: 


k? _ k 4 

72 » 74 - — — 

a " 2 


7 4 = 74"n 


•q = ( 1 -e" 2 ) l/2 ; 0 = cosi" 


A = 9 : A, = 1 (1 - 110 2 -4O0 2 A 0 ) 

° T~ 1 ^ 

(1 -50 2 ) 8 


2 + e " 2 - 400 e" 2 0 2 A 2 0 


A 2 = Ti 3 7' 2 A,-_i.7'2 < - 40(5e" 2 + 2)0 2 A 0 


16 


110 2 (3e + 2) 


74 = ; 7'j = 73^' 

‘ «" 3 


A, = - -L 0(H + 80 A 0 + 200 A 2 ) 
8 


A 4 = 


i ( 1± ) sini"; A 5 = e"0 4 _ 

4 \ 7 r 2 / 0 + 0) 


= n„ < 


7 2"n 


A (30 2 - 1) + ~l ' 2 | iL ] [0 2 (3O-96 ti-9Oti 2 ) 


+ (16x| + 25-n 2 - 15) + 0 4 (144 ti + 25j\ 2 + 105)] 


+ ( i£ ) Tie" 2 7' 4 (3 + 350 4 - 3O0 2 ) 
16 


n., < 


3 7S(50 2 -1)+ A 7' 2 [25'0 2 + 24ti-35 
2 32 2 

+ (90- 192-n- 126-n 2 )0 2 + (385 + 360 t) + 45ti 2 )0 4 ] 

+ A 7' 4 [21 -9-n 2 + (126 t! 2 -27O)0 2 + (385- 189t] 2 )0 4 ] 
16 


h" = n 0 < 


- 37 2 9 + 

A 7' 2 [(9tT + 1 2tj - 5) + (- 5-r| 2 - 36 t| - 35)0 3 1 


+ i_7'4(5-3Tf)0(3-70 2 ) 


2 



Compute secular terms: 


V* = “mean'’ mean anomaly 
r = r o + (n 0 + t") (t - 1 0 ) 
g" = mean argument of perigee 

g" = g"o + g" (t - to) 

h" = mean longitude of the ascending node 

h" = h" 0 + h" (t - t„) 

u" = “mean" mean argument of latitude 

u" = r + g " 

Compute long period terms: 

5,e' = r| 2 A 4 sin g" + e" t| 2 y' 2 A, cos (2 g") 

5 2 e' = e" T| 3 y' 2 A, sin (2 g") - -q 3 A 4 cos g" 

e' = y^e') 2 + (e" + 8 ie ') 2 

x • -e"8,e' 

o,i = — 

Ti 2 tani" 

h' = h" + e " 2 A, y 2 sin(2g") + (e"0 —S- )cosg" 

sin~i" 


u' = u" + A 2 sin (2 g") + [(2 + iq - e" 2 )e" + A 5 ] cos g" 

(l+iq) 

t'=tan l 8 2 e'cos/" + (e"+8,e')sin/" 

_ (e"+ 8,e')cos/"- (8 2 e'sin/") 

g' = u' - 


however for e"< 0.05 we define: 

g' - g" and (' = t" 


and continue with, 


E' - e' sin E' = /' 



f ' = tan 


.1 [ V l-e' 2 sinE' \ a" 


cos E' - e' 


’ a - 17 
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r' 1-e'cosE' 


Semi-major Axis 


a = a" {1 -h -y 2 [(30 2 - 1) (a 3 - t) _3 ) + 3(1 -0 2 ) a 3 cos (2g'+2f)]} 


Compute eccentricity: 


N .1-2 


f 3 — 7 2 (1 - e 2 )cos(2g' + 2f')[3e"cos 2 f' + 3 cos f' + e" 2 cos 3 f' + e"] 


8,6 = 8^' + — if < 


- y f 2 (1 - 0 2 ) [3 cos(2g' + f') + cos(3f' 4- 2g')] + 


. ,3e: - ,hj (^) [ ^ + ( T ^) 


4- 3e" cos" f 4- 3 cos f 4- e 1 


n2 


S^e^Soe'- -ti 3 7' 2 
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f 2(30 2 - 1) (a 2 ti 2 + a + 1) sin f 

(-a 2 r| 2 -a+ l)sin(2g' + f') 


+ 3(1 -0 Z ) 


+ (a 2 Ti 2 + a + I ) sin (3f' + 2g') 
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\/ (82c)" + (e" + 8,e)“ 


Compute inclination: 


i = i" + 8,i+ I 7' 2 0 V 1 -0 2 [3 cos (2 g' + 2 f') + 3 e" cos (2 g' + f') 


+ e" cos (2 g' + 3 f')] 


Compute longitude of the ascending node: 

o 1 / 


h = h' - - 7 ' 2 0[6(f' -l' + e" sin f') -3 sin (2g' + 2f ) - 3 e" sin(2 g' + f') 
2 


e" sin(2g' + 3f')] 


> 

cos 3 f'] 
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Compute mean argument of latitude. 



r 

( i + a 2 ti 2 + a)sin (3f' + 2g') 

u-u'+-( \ e"y' 2 in . 

4\ti+W 

3(1 - 0 2 ) 

1 

3 

_ + ( 1 - Tt 2 a 2 - a) sin (2g' + f') 


[ + 2(30 2 - 1 ) (T| 2 a 2 + a + 1) sin f 


+ 


+ 


--y' 2 (5 0 2 - l)(e" sinf' + f -/') 

2 

1 2 f e"sin(2g' + 3f') + 

_ y 2 (3 - 50 ) ^ 3 (2 g , + 2 f ) + e " sin (2g' + f')] 
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Compute mean anomaly: 


t =tan -1 


8 2 e cos l " + (e" + 8]e) sin l " 
(e" + 8 , e) cos / " - (8 2 e sin 1 ") 


Compute argument of perigee: 


g = « - 1 


Stable Osculating to Mean Conversion 

Walter’s 11 algorithm for osculating to mean conversion is unstable for low e in Keplerian space; the apparent 
instability of the § iterative osculating to mean element conversion is removed by translating the iteration from 
mean Keplerian space to mean Cartesian space. 


Define: 

Q = ( a " ? e ", i", g", h", l ") - Mean Keplerian Elements 
0 = (a,e,i,g,h,<f) - Osculating Keplerian Elements 
X = ( x ", y", z", x", y", z") - Mean Cartesian State 
Elements 

Y = ( x , y, z, x, y, z)- Osculating Cartesian State 
Elements 


Given an osculating Cartesian state Y we determine 

0 (o) — f 2B (Y) 
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Where f 2B represents the Keplerian state two-body functional 
the iterative algorithm. 


relationship to the Cartesian state. Then employing 


£ (J) — B g (0 <j) , At = o) 

Y 0> — f^ B (0 a> ) 

X i (j+1) = Xf + (Y s - Y ( f), i = l, 2 6 

n (j+i) — f 2B (x (j+i) ) 


For j 0, 1,2, ..., 10 or until the following criterion is satisfied: 


I Y, - Y<j> 


Where e is some preassigned small 
relationship. 


positive number. Let this algorithm be represented by the symbolic functional 


Q *— O(Y) 


Semimajor Axis Decay Rate 

deca^M^““"fw“h" ,ean C0 " VerSi ° nS “ mea " Peri ° d ' F) ,nMrValS ' We de '™ e ,he Sem,ma J OT “ is 


■* Oj(Y) 


given for i - 1,2 M; we compute the mean semimajor axis decay rate as: 


a = 



Orbit Propagation With “Bg” 

To update the orbital elements to time (At = t - t G ) with the Bg-theory, we assume the orbital elements re- 

<P) and rec,lf> the lheo,, '' s cons,ants a< one period imervais ,wi,h ,he 5 **** ra «> 


N = int 
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Employing the above iterative method, we have: 


a j = 


j-i 


+ aP 


j-i 


with perigee constant, the eccentricity decay equation is given by. 


e" = e"., + _ — —L1L- aP 




a j-i 


■H 


(O = (Oj-i - 7 f ) aP 2 .,^ 


4 V a: 


j-i 


GM 

q"3 


Evaluating only the secular part of Bg, we obtain: 

n j *-B"g(a j ,p j . l) 


Where j =1,2, N = > 


On (a n , e n , j" N , g" N , h" N , /" N ) 


then at time At the osculating elements are given by evaluating the full Bg theory with: 

ft (At) — Bg (ft(T N ), At - T n ) 

With T n = t c + N x P. Let us represent the semianalytic theory with the rectification algorithm for retarded mo- 
tion symbolically by Bg. 


0(J 2 ) Secular Retarded State Matrizant 

Spacecraft state estimation algorithms employing either a weighted least square estimator or sequential filter 
estimator for orbit refinement and prediction requires development of a set of partial derivatives called the matri- 
zant, or state transition matrix. These partial derivatives give the relationships between perturbations in the 
spacecraft state at observation times to perturbations in the state at epoch. 

The objective of this study is to obtain an analytical formulation to employ with the Bg theory which nre- 
serves its accuracy and computational efficiency for micro-computer application. As Rice 12 has noted, neglecting 
small forces such as gravitational harmonics and drag may result in non-negligible errors, and “a good approxi- 
mation to the state transition may be obtained with a simplified force model provided that the transition matrix is 

tr r aJeC u? ry baSCd ° n the S3me force model ” However, for retarded artificial satellites orbiting 
near the Earth, the J 2 oblateness potential and semi-major axis decay rate contributes strong perturbations which 
may not be neglected for accurate satellite orbit predictions. 

With the foregoing observations in mind we differentiate the secular solution of “Bg” with J, only to obtain 
t e transition matrix in terms of the keplerian elements and the empirically determined semi-major axis decay rate 

parameter (a). If the state (secular retarded keplerian elements) is given as a function of the initial state ft and 
time(t)i.e., — ° 
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0(t) = ft(ft 0 (t 0 ), t) 


then <f> (the matrizant) may be determined by taking partial derivatives of this functional relation: 

<KM„) - ^ (Qo(to).t) 


with an element of <P(t, t Q ) given by 


4>ij(t, U 


dQj(t) 

dOj(U 


where (i, j = 1, 2, .... 7); ft is a generic symbol for the dynamic orbital elements and the decay rate parameter of 
a state vector of dimension 7. 


Computation of Matrizant Elements 


Define 


| J 2 R 2 ;P = a(l -e 2 );n 



with: 

a = a„ + a At 
e = e„ 4- ^ ~ e - a At 

a o 

' = io 

g = go + g^t 

h = h G + hAt 

7 = 1 0 + nAt+ 1 2 At- - — a At 2 
4 a 0 

g = - 7 (5 cos 2 i 0 - l)n 

P 

h = -2 cos i 0 n 
P 

7 2 = (3 cos 2 i Q - 1 )n 
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then the partials are: 
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Simulated Trajectory Data 

Simulated trajectory data, i.e., osculating state vectors were used to simply demonstrate the Bg-theory capa- 
bility to represent retarded satellite motion about an oblate planet. This is not meant to impute any claim of accu- 
racy to the method. Fig. 1 presents a “Truth Ephemeris” generated by the Cowell numerical propagation of near 
circular 400 km altitude satellite perturbed by a 4 x 4 gravitational field and drag, compared with the ephemeris 
generated by the Bg-theory without drag (i.e., the retardation rate parameter is set to zero). An error growth of 
1000 km in along-track and 80 km in radial is realized over a 6-day span. The constants of the Bg-theory (with 
the decay rate included) is derived from a fit over a 3-day span to the “Truth Ephemeris” state vectors by a diffe- 
rential correction technique. The post fit error growths are reduced to less than 2 km (Fig. 2). 

Tracking Data 

Real tracking data demonstrates the Bg-theory is favorable orbit determination and prediction capabilities. 
Orbit determinations for a number of different epochs employ real (SMM-Solar Maximum Mission) tracking data 
over a 2-day span to differentially correct the epoch state and drag model constant for the Cowell method and the 
epoch mean elements for the Bg-theory. The predicted ephemeris of the Cowell method and the Bg-theory is then 
compared with a series of definitive state solutions determined over successive 2-day spans. The comparable re- 
sponses of the Cowell and Bg methods are presented in the table. The table demonstrates that for all practical pur- 
poses, wherein drag is significant, the theory is in fact as accurate or non-accurate as a Cowell method in a post 
differential-corrected orbit ephemeris prediction. 
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a" = 6775.8813 
e” = .00057510273 
i" = 28.78258 
g = 247.91408 
h" = 19.780076 
C"= 112.08581 
DOT a - 0 


MAXIMUM-km. 

ALONG = 1067.99 
RADIAL = 82.801 
CROSS = 0.636 
TOTAL = 1071.195 


a"= 6775.9246 
e"= .00062668228 
i" = 28.782232 
g" = 243.5946 
h" = 19.778049 
116.4108 

DOT a = -4.71 27509E-6 


MAXIMUM-km. 

ALONG =- 1,906 
RADIAL = 1.021 
CROSS = - 0.302 
TOTAL = 1.993 


ALONG track 



ALONG TRACK 



RADIAL 



Fig. 1 Cowell v. “Bg” (a =0) 



t 1 2 3 4 5 6 7 8 days 

6 hr 


Fig. 2 Cowell v. “fig” (a * 0) 
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Table 1 . SMM In-Track Errors (km.) 



COWELL 

Bg 

COWELL 

Bg 

COWELL 

Bg 

EPOCH 

80/o: 

1/14 

80/03/20 

80/08/1 




D.C. 

ARC 




1 

mm 

■9 


0.09 

0.71 

2 


warn 

u 


0.12 

0.85 


PREDICTS 




2 

-0.44 

0.66 

-2.96 

-1.87 

2.41 

2.31 

4 

-2.84 

-0.81 

-6.78 

-4.08 

6.43 

4.67 

6 

-10.39 

-7.46 

-19.62 

-14.66 

13.02 

8.73 

8 

-27.15 

-23.28 

-42.60 

-34.60 

22.88 

15.28 


-51.66 

-46.88 

- 77.20 

-65.26 

38.18 

26.60 

12 

-91.94 

-86.26 

-123.75 

-106.96 

62.04 

45.88 

14 

-149.09 

-142.39 

-176.93 

-154.47 

95.58 

74.00 

16 

-224.58 

-216.62 

-237.02 

-208.23 

141.93 

114.11 


Conclusion 

A computationally efficient semianalytic orbit theory for satellite motion perturbed by oblateness and drag ef- 
fects with minimal demand on computer storage requirements has been developed. A stable osculating to mean 
conversion algorithm is presented which is used to provide an accurate first order estimate to the semi-major axis 
decay rate for Bg. This algorithm can also be used to produce a reference orbit of mean trajectory state vector 
data for a stable adaptation of the constants of the theory to accurate special perturbation methods, e.g., Cowell, 
by differential correction. Also a computationally efficient 0(J 2 ) secular retarded state matnzant is developed to 
employ with the theory in batch or sequential estimation algorithms for orbit refinement. The paper presents com- 
plete simplified computational algorithms for ephemeris generation and state estimation. 

Figures 1 and 2 graphically demonstrate how well the semi-major axis decay rate along with the period ratifi- 
cation algorithm accommodates significant drag effects. The table establishes that for all practical purposes t e 
theory is in fact as accurate or non-accurate as a “definitive Cowell method” in post differential-corrected orbit 

P Furthermore, suggestions are made for the inclusion of J 3 long period computations that can significantly im- 
prove the accuracy of existing programs of the Brouwer or Brouwer- Lyddane theory, and for resolving the radial 
error discrepancies between the respective theories. 

The theory has been implemented on a very inexpensive personal home micro-computer with accuracy com- 
patible to that obtained on the large mainframe IBM 360-95 computer. The power supply for the T/S 2068 re- 
quires only 1 ampere. Given the computational efficiency, low storage requirements, and low power consump- 
tion, such economized general-perturbation models could be utilized as backup orbit propagators onboard a fu- 
ture spacecraft. The theory’s practical use for portable field applications at the foot of an antenna has been de- 
monstrated, e.g. , tracking a satellite to control station pass acquisitions. ^ 

Furthermore, the economized theory lends itself to iterative closed analytic algorithms (Gordon ) for au- 
tonomous onboard future event predictions; e.g., occultations, station contacts, special ground trace reference 
points, etc. This eliminates costly numerical techniques which require superfluous evaluations of constraint 
equations. 
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A V* ec ° nomized general-perturbation models discussed in the paper retain the physics of the retarded 
P rob| em and require a set of only seven constants for initialization. This provides superior prediction accu 
a differed expen encing significant drag and oblateness perturbation effects which de^ay gracefully beyond 

arc - w,,h cons,derabie da,a >— • — — .ss 

.ion E S yaems sucks' OTDFS m 7™'? supporttd ** me " s ' K «etem,i M - 

r ra G m \rT ssca r co,,a ^^ 

fccur^tr eph7mS ^ COmpura,lo " al «*<*«» and sufflcL, 

This should facilitate the theory’s use: for alternative onboard orbit propagation or as a back-up propagator 
with mobile or fixed ground tracking installations; by experimenters; by satellite project scientists- f e a^y user 
ho has need of ephemens data readily accessible, cheaply obtained, and where routinely generated. ’ Y 
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Appendix 

-'Bt-Mhfintv “ S,ing “ TS 2 a°°° Basic l«0) for satellite ephemens genera, ion w,,h the 

j. J rX , - °fo |ls generated on the T/S 2068 or the large main frame IBM 360-95 computer yielded ran 

& y me rfg theory on the T/S 2068 is presented Appendix R nrp^nts a RA^tr^ ♦ . i 

dtffercmial (correction ,DC, program utilizing the P 'Bg” theor/S ^ ^ 

sen ed in the paper tor spacecraft state estimation by processing trajectory data as input from a reference orhif a 
sample ouipu, of a differentially connected orbit is prcLted which 

s. However, degradation in the DC is more pronounced on the T/S 2068 PC than using double precision 
anthmettc on , he large main frame IBM 360-95. As noied, ihete is no apprcciable difference hf 0^4^ " 
ation. Appendix C presents a Basic listing and sample output of the stable osculating to mean elements conver 
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Appendix A 




mi! SS V i;Sffici; s..!-*«ly«i=.l orbit Th.or, 

r» nis - «" I rnsr 1 *’ " D - 

AIAA 24th Aerospace Sciences Meeting Jan. 6-9, 1986 

T * 

5 DEF FN M(X,N)=X-N*INT (X/N) 

8 SS fS ^C.^S^ANd’siO^-I AND S<0 OR ♦ . AND S>=0>,»CS 


9 CLS 

10 GO SUB 2900 
20 REM EPHEM 

22 FOR t=0 TO ENDt STEP DELt 

23 LET nt=t/DELt+l 

25 GO SUB 270 

26 GO SUB 300 

27 GO SUB 80 
29 GO SUB 90 

37 NEXT t 

38 BEEP .5,0: BEEP 1,12 

39 STOP 

80 REM t-Day ,hr , min, sec 

81 LET Day = I NT Ct/86400) 

82 LET hr = I NT ( < t-Day*86 400 ) / 3600) 

83 LET m i n= I NT < < t- < Day*86400+hr*3600 ) > /60> u 

84 LET sec= I NT ( t- C Day*86400+hr*3600+m l n*60 ) ) . LET 

c+ . 5 > 

86 PRINT INVERSE 1 ; Day 5 " Day _ " i hr 5 ; hr " ; 5 seCJ ” S 

87 LPRINT Da y ; ” Da y “ : hr ; hr ,min, min ,s , 

88 PRINT 

89 RETURN „ _ _ 

90 REM PRINT a , e , i , g , h , 1 -x , y , z , Dx , Dy , Dz 

91 PRINT AT 0,20;” ^-km/sec ” ' [ . , . . . . t * RTD , - z = 

92 PRINT a= la, x= ,x, km •-»*»?;' - D z=”;Dz 

.-i 9 .RTD,-DK--!Dx--h--lh*RTD,-Dy.-iDy 1-_.1*RTD, Dz ,m 

93 LPRINT l 6 ‘ "I:: 1 1’ “* 6 . “h" - ;JrtdI 

i *RTD ; TAB 16; “z = » z g= , g*RTD , TAB 16, Ox 

; -Dy=“ ;Dy’ "1=" ; 1 *RTD ; TAB 16; “Dz=” ;Dz 
99 RETURN 

100 REM SOLVE KEPLERS EQ . 

1 10 LET EA=0 

115 IF 1=0 THEN GO TO 160 

‘,25 PM SiVro 10: LET OEA-EA: LET FE-EA-e.SIN EA-1: LET 
FE/(l-e*COS C EA-0 . 5*FE) ) : LET DEA = ABS (EA-OEA) 

135 IF DEA<=0 . 1 E-8 THEN GO TO 160 
140 NEXT N 

160 LET EA=FN M(EA,2*PI> 

199 RETURN 

20? S aS?°"“ T LET EDP.,0: LET IDP-iO. LET =DP= 9 0: LET H 
LET LDP= 1 0 

203 L K EDP2-EDP*EDP? P LET^CN2=N 1 -EDP2 : LET CN=SQR <CN2> 

204 LET GM2 = K2 / ADP~N2 : LET GMP2 = ^' < £N2*CN2 > : LET GM4- < 
4 : LET GMP4=GM4/CN / '8 : LET F1MG2-F D4*GM p 2 

205 IF Dt = 0 THEN LET CI=COS (IDP): LET CI2-CI*L1. 

Cl : LET Cl 4=CI2*CI 2 

206 REM 1D0T , gDOT , hDOT 


" e= “;e,"y=";y’“i = "; i*RTD, “z= 
3y = ” ; Dy ’ " 1 = “ ; 1 *RTD, "Dz=” ;Dz 
;x;“km“'"e=";e;TAB 16;“y="n 
: TAR 16; "Dx = “ ;Dx’ "h=* :h*RTDi 
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207 LET 1D0T=N0*(CN*(GMP2*(F3D2*(N3*CI2-N1)+F3D32*GMP2*(N25*CN2 
+16*CN-15+(N30-96*CN-N90*CN2)*CI2+(N25*CN2+144*CN+105)*CI4)))+F1 

5D 1 6 *GMP4*CN*EDP2 * (N3-N30*CI 2+N35*CI 4) ) 

208 LET nD0T=N0+ 1 DOT : LET PD=PI2/nD0T 

209 LET gDOT=NO*((GMP2*(F3D2*(N5*CI 2-N 1 ) +F3D32*GMP2*(N25*CN2+24 
*CN-N35+(N90- 1 92*CN- 1 26*CN2)*CI 2 + ( 45*CN2+360*CN+N385)*CI4) ) ) +F5D 
1 6*GMP4* ( 2 1 -N9*CN2 + (N 1 26*CN2- 270 ) *CI 2+ ( N385- 1 89*CN2 ) *CI 4) ) 

210 LET hDOT = NO*< < GMP2* ( F3D8*GMP2*( (N9*CN2+ 1 2*CN-N5)*CI - ( N5*CN2 
+36*CN+N35)*CI 3)-N3*CI)) +F5D4*GMP4*(N5-N3*CN2 ) *CI*(N3-7*CI2)) 

211 RETURN 

220 REM SECULAR-GDP , HDP , LDP 

221 LET GDP=gO+gDOT*Dt 

222 LET GDP=FN A(GDP) 

223 LET HDP=hO+hDQT *Dt 

224 LET HDP=FN A(HDP) 

225 LET LDP=10+nD0T*Dt 

226 LET LDP=FN A (LDP) 

227 LET a = ADP : LET e = EDP : LET i = IDP: LET g=GDP : LET h=HDP : LET 
1 =LDP 

229 RETURN 

230 REM SP, LP- CONSTANTS 

232 LET CN3=CN2*CN : LET CN6=CN3*CN3: LET F ID 1 CN= 1 / ( 1 +CN) : LET F 
IDCN3=1/CN3: LET F1DCN6=1/CN6 

233 LET GM3=K3 / ADP" 3 : LET GMP3=GM3/CN6 : LET G3DG2=GMP3/GMP2 

234 IF Dt=0 THEN LET SI=SIN (IDP): LET TI=SI/CI: LET P3T2M1=N3 
*CI 2-N 1 : LET P 1 MT2 = N 1 -Cl 2 : LET SQ1MT2=SQR (P1MT2): LET T31MT2 = N3 
*P 1 MT2 : LET T5T2M1=N5*CI2-N1 : LET P3M5T2=N3-N5*CI 2 : LET A0=CI2/( 
N 1 -N5*CI 2 ) : LET A 1 =F 1 D2*F 1 D4* ( N 1 -N 1 1 *CI 2-N40*CI 2*A0) : LET A3=-F1 
D2*F1D4*CT*(N1 1 +80* A0+ 200* AO* AO) 

235 LET EDPT3=N3*EDP : LET SP3=F 1 D2*GMP2 : LET TSP3=CI*SP3: LET S 
P6=CI *SP3*SQ 1 MT2 

236 LET A2=CN3*GMP2*A 1-F1D4*F1D4G2*( N2+EDP2- 400*EDP2*CI 2*A0*A0- 
40*(N5*EDP2+N2 )*CI 2*A0- 1 1 *C I 2*(N3*EDP2+N2) ) : LET A4=F 1 D4*G3DG2*S 
I: LET A5=(A4*EDP*CI ) /(Nl+CI ) 

239 RETURN 

240 REM UDP, PERIODIC TERMS 

241 LET EP=EDP : LET GP=GDP: LET LP=LDP: LET UDP=GDP+LDP: LET UD 
P=FN A(UDP) 

242 REM LP-TERMS 

243 LET SG=SI N (GDP): LET CG=COS (GDP): LET S2G=N2*SG*CG : LET C 
2G=N2*CG*CG-N 1 

244 LET D 1 E=A4*SG+EDP*GMP2*A 1 *C2G : LET D 1 1 =- ( EDP*D IE) /TI : LET D 
1 E=CN2*D 1 E : LET D2E=EDP*CN3*GMP2*A 1 *S2G-CN3*A4*CG 

245 LET EP=SQR (D2E*D2E+ (EDP+D 1 E ) *( EDP+D IE)) 

246 LET HP=HDP +EDP2*A 3*GMP2 *S2G+ ( ( EDP*C I*A4)/(SI*SI))*CG: LET H 
P=FN A ( HP ) 

247 LET UP=UDP+A2*S2G+ ( ( EDP*A4*F 1 D 1 CN ) * ( N2+CN-EDP2 ) +A5 ) *CG : LET 
UP=FN A ( UP ) 

248 LET SL=SI N (LDP): LET CL=COS (LDP) 

249 IF EDP > = 0.05 THEN LET SM = D2E*CL + ( EDP+D 1 E ) *SL : LET CM=(EDP + 
DIE) *CL-(D2E*SL) : IF CMOO THEN LET TM=SM/CM: LET LP=FN 1(SM,CM 
• TM ) : : LET GP=UP-LP: LET GP=FN A(GP): LET SG=SIN (GP): LET CG=CO 
S (GP): LET S2G=N2*SG*CG : LET C2G=CG*CG-N1 

250 REM FP 

251 LET 1 =LP : LET e=EP: GO SUB 100: LET EAP=EA : LET SEA=SIN (EA 
): LET CEA=COS (EA) 

252 LET ADR=N 1/(N1-EP*CEA): LET ADR2=ADR*ADR : LET ADR3=ADR2*ADR 
: LET SF= ADR*SQR ( N 1 -EP*EP ) *SEA : LET CF= ADR* ( CEA-EP ) : LET FP=FN 
T(SF, CF) 

253 REM SP-TERMS 

254 LET CF2=CF*CF: LET CF3=CF2*CF: LET S2F=N2*SF*CF: LET C2F=N2 
*CF2-N1: LET S3F=N3*SF-N4*SF*SF*SF : LET C3F=N4*CF3-N3*CF : LET S2 
GPF=S2G*CF+C2G*SF : LET S2GP2F=S2G*C2F+C2G*S2F : LET S2GP3F=S2G*C3 
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F+C2G*S3F: LET C2GPF=C2G*CF-S2G*SF : LET C2GP2F=C2G*C2F-S2G*S2F : 
LET C2GP3F=C2G*C3F-S2G*S3F 

255 REM COMPUTE a,e,i,g,h,l 

256 LET a= ADP* (N1+GM2*(P3T2M1 *( ADR3-F 1DCN3) +T3 1MT2*ADR3*C2GP2F ) 

) 

257 LET D1E=(F1D2*CN2*( < N3*F 1 DCN6*GM2*P 1 MT2*C2GP2F*(EDPT3*CF2+N 
3*CF+EDP2*CF3+EDP ) ) - (GMP2*P 1 MT2*( N3*C2GPF+C2GP3F) ) +P3T2M 1 *GM2*F 1 
DCN6*(EDP*CN+EDP*F1D1CN+EDPT3*CF2+N3*CF+EDP2*CF3> ) )+DlE: LET D2E 
=-FlD4G2*CN3*(N2*P3T2Ml*(ADR2*CN2+ADR+Nl ) *SF+T3 1MT2*< < -ADR2*CN2- 
ADR+N l ) *S2GPF+( ADR2*CN2+ADR+F 1 D3) *S2GP3F ) )+D2E: LET e =SQR <D2E*D 
2E+C EDP+D 1E)*(EDP+D1E) ) 

258 LET i=IDP+DlI+SP6*(N3*C2GP2F+EDPT3*C2GPF+EDP*C2GP3F) : LET i 
=FN A ( i ) 

259 LET h=HP-TSP3*( N6* ( FP-LP+EDP*SF ) - ( N3*S2GP2F +EDPT3*S2GPF +EDP 
*S2GP3F ) ) : LET h = FN A(h) 

260 LET u=UP+ ( F 1 D 1 CN*F 1 D4G2*EDP*CN2*( T3 1 MT2*< S2GP3F*< F 1 D 3+ ADR 2* 
CN2+ADR)+S2GPF*(N 1 - ( ADR2*CN2 + ADR > > ) +N2*SF*P3T2M 1 *< ADR2*CN2+ADR+N 
1 ) ) )+GMP2*F3D2*(T5T2Ml*(EDP*SF+FP-LP> >+P3M5T2*(FlD4G2*(EDP*S2GP3 
F+N3*(S2GP2F+EDP*S2GPF) ) ) : LET u = FN A(u) 

261 LET SM=D2E*CL+ ( EDP+D IE) *SL : LET CM=(EDP+D1E)*CL-D2E*SL: IF 
CMOO THEN LET TM = SM/CM 

262 LET 1 =FN 1(SM,CM,TM) 

264 LET g=u- 1 : LET g=FN A(g) 

269 RETURN 

270 REM ORBGEN 

271 LET Dt=t : IF Dt=0 THEN GO SUB 200: GO SUB 230 

272 IF D0Ta=0 THEN GO TO 290 

273 IF t=0 THEN LET tsum=0: LET t0=0: GO TO 289 

274 LET s i gn= 1 : LET Dt = t-tO: IF Dt<0 THEN LET siga = -l 

275 LET PD=sign*PD 

276 LET nPD= 1 : IF ABS Dt > = ABS PD THEN LET nPD=INT <Dt/PD) 

277 IF ABS t < ABS (tsum + PD) THEN GO TO 289 

278 FOR n= 1 TO nPD 

279 LET DOTe = ((l-eO)/aO) *DOTa 

280 LET D0TnD2= - C 3 / 4) * ( nO/aO) *DOTa 

281 LET DELa=DOTa*PD : LET DELe =DOTe *PD : LET DELI =D0TnD2*PD*PD 

282 LET aO=aO+DELa : LET eO=eO+DELe: LET 10=I0+DELI: LET 10=FN M 
(10, 2*PI ) 

283 LET Dt=PD: LET tsum=tsum+PD 

284 GO SUB 200: GO SUB 220 

285 LET aO=a : LET e0=e: LET gO=g: LET hO=h: LET 10=1 

287 NEXT n 

288 GO SUB 200: GO SUB 230 

289 LET tO=tsum : LET Dt=t-tO 

290 GO SUB 220: GO SUB 240 

299 RETURN 

300 REM KEP-POSVEL 
302 GO SUB 100 

304 LET cosEA=COS EA : LET sinEA=SIN EA 
308 LET e 1 =a*SQR (l-e*e> 

310 LET r=a*( l-e*cosEA) 

312 LET DotEA=SQR (GM/a)/r 
314 LET Xw=a* ( cosEA-e ) 

316 LET Yw=el*sinEA 
318 LET DXw=-a*DotEA*s i nEA 
320 LET DYw=e 1 *DotEA*cosEA 
322 LET sini=SIN i: LET cosi=COS i 

324 LET s i ng=S I N g: LET cosg=COS g 

326 LET s inh = SIN h: LET cosh = COS h 

328 LET Px=cosg*cosh-s iag*s inh*cos i 
330 LET Py =cosg*s i nh+s i ng*cosh*cos i 
332 LET Pz=sing*sini 

334 LET Qx=-s i ng*cosh-cosg*s i nh*cos i 
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336 LET Qy = - s i ng *s inh+cosg*cosh*cos i 

338 LET Qz=cosg*sini 

339 REM x,y,z 

340 LET x=Xw*Px+Yw*Qx 
342 LET y=Xw*Py+Yw*Qy 

344 LET z=Xw*Pz+Yw*Qz 

345 REM Dx,Dy,Dz 

346 LET Dx=DXw*Px+DYw*Qx 
348 LET Dy=DXw*Py+DYw*Qy 
350 LET Dz=DXw*Pz+DYw*Qz 
399 RETURN 

2900 PRINT " INPUT-BgORBIT" 

2901 GO SUB 9900 
2905 PRINT “INPUT aO“ 

2910 PRINT "INPUT eO" 

2915 PRINT "INPUT iO" 

2920 PRINT “INPUT gO" 

2925 PRINT “INPUT hO“ 


INPUT 

al : 

PRINT 

al 



INPUT 

e I : 

PRINT 

e I 



INPUT 

il : 

PRINT 

i I : 

LET 

i I = i I *DTR 

INPUT 

gl : 

PRINT 

gi : 

LET 

g I =g I *DTR 

INPUT 

hi : 

PRINT 

hi : 

LET 

hi =hl *DTR 

INPUT 

1 I : 

PRINT 

1 I : 

LET 

11=11 *DTR 


INPUT DOTa: PRINT DOTa 


2930 PRINT “INPUT 10" : 

2935 PRINT "INPUT DOTa 

2946 PRINT "INPUT OUTPUT DEL ( t ) in Sec’s 
2948 INPUT DELt: PRINT DELt ; “Sec ’ s “ 

2950 PRINT "INPUT OUTPUT SPAN IN hrs 
rs": LET END t= ENDt *36 00 
2970 LET a0=al : LET e0=el: LET i0=il 

10=11 

2975 LET t=0: LET tsum=0 
2977 PRINT "a 
" ;hI*RTD' "l ’ ’ 


INPUT ENDt: PRINT ENDt ; " h 
LET g0=gl : LET h0=hl : LET 


■ ’ ’ = " ; a I ’ " e ’ ' = “ ; e I ’ " i ’ ’ : 
= " ; 1 1 *RTD ’ " DOTa= " ; DOTa 


2978 LPRINT "a 


» _ f * — 


“ ;al ’ “e ” =“ lei 


i 


“ ; i I *RTD ' "g’ ' = “ ; gl *RTD’ “ h ’ ' = 
= "; iI*RTD’“g’ ’ = “ ; gl *RTD’ ” h ’ ' 


= ; hi *RTD' “ 1 ’ ’ = “ ; 1 1 *RTD’ “ DOTa= " ; DOTa 
2990 PRINT “RETURN": PAUSE 300 
2999 RETURN 
9900 REM "CONSTANTS" 

•9901 RESTORE 

9920 DATA 39860CK6 3^ 6 378 . 1 66 ! 0 . 7 29 2 1 1 59E- 4 , 298 . 25 , -0 . 1 08265 1 7E-2 

492rREArF?D2:FiD3!FlD4!F3D2,F3D8,F3D32,F5D4,F5D16,F15D16 

9926 READ & N^S ! Si 1 f I SSSTSai I SSoI Si2?i iS . -BBS 

9928 DATA 1,2,3,4,5,6,9,11,25,30,35,40,90,126,385 
9930 LET Ke =SQR (GM/Re^S) 

9940 LET K 2 = -0.5*J2*Re / '2 
9950 LET K3= J3*Re ^ 3 
9960 LET K 4 = F3D8*J4*Re / '4 
9970 LET P I 2= 2*PI 

9980 LET RTD=180/PI: LET DTR=PI/180 
9999 RETURN 
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a’ ’ =6775 .8813 
e’ ’=.00057510273 
i ' ’ =28.78258 
g' ’ =247 . 9 1 408 
h’ ’=19.780076 
1 ’ ’ = 1 12.08581 
D0Ta=0 

ODay Ohr Omin Osec 
a=6778 .1401 x=637 1 . 8445km 

e= .0010000099 y=2291.4926 

i = 28 . 8 z = - . 00017868567 

g= .00012070225 Dx=-2 . 2763895 

h= 1 9 . 78 Dy=6. 3298482 

1=359.99988 Dz=3. 6980524 

ODay 2hr Omin Osec 
a = 677 4 . 067 3 x =- 386 1 . 6402km 

e=. 00078805481 y=4619.9292 

i=28.768669 z=3092.0373 

g=117. 97021 Dx = -6 . 1 8656 1 1 

h=19. 166926 Dy = - 4 . 39 390 1 8 

1=350.39302 Dz=- 1 . 1 6 35038 

ODay 4hr Omin Osec 
a=6776 . 53 x = - 40 1 0 . 760 1 km 

e=. 001 1819574 y=-5092.5391 

i=28. 787583 z=-1947.4903 

g=216. 38997 Dx=6. 0713698 

h=18. 631543 Dy=-3 . 647801 3 

1=0.29937848 Dz = - 2 . 965 1 278 

ODay 6hr Omin Osec 
a = 6776 .6611 x = 6329 . 5807km 

e=. 00091255666 y=-1515.9705 

i=28. 788601 z=-1865.5031 

g=311. 04825 Dx=2. 4767267 

h=17. 971268 Dy=6. 6030454 

1=14.027272 Dz=3. 0314147 

ODay 8hr Omin Osec 
a=6773 . 988 1 x= 1 44 . 1 8588km 

e= .00068098079 y=6005.1244 

i=28. 768061 z=3121.8257 

g=93. 110521 Dx=-7 . 5938972 

h=17. 433423 Dy=-0 . 3686 32 47 

1=340.2922 Dz=l. 056009 

ODay lOhr Omin Osec 
a = 6778 . 1 385 x = -6423 . 1 55 lkm 

e=. 001387741 1 y=-2132.785 

i=28. 799978 z = -100. 26728 

g=182. 50246 Dx=2. 1729597 

h=16. 824348 Dy=-6 . 3698705 

1=359.26165 Dz=-3 . 6977325 

ODay 1 2hr Omin Osec 
a=6774 .1573 x=3794 . 6973km 

e= .0007266 17 1 4 y=-4698.0551 

i=28. 769325 z=-3058.5189 

g=271. 64378 Dx=6. 2734922 

h=16. 207078 Dy=4. 2364246 

1=18.485857 Dz=l. 2722066 

ODay 1 4hr Omin Osec 
a=6776 . 3894 x = 4 1 09 . 8 495km 

e= .0008874069 y=4984.1467 

i=28. 786536 z=2026.628 

g=50. 552286 Dx=-6 . 0 1 87335 

h=15. 674014 Dy=3. 7838106 

1=347.90236 Dz=2. 8951282 


a’ ’=6775.9246 
e ” = . 00062668228 
i’ ’=28.782232 
g’ ’ =243.5946 
h’ ’ =19 .778049 
1 ’ ’ = 1 16 . 4 108 
D0Ta=-4 . 7 1 27509E-6 
ODay Ohr Omin Osec 
a=6778 . 1 825 x=6372 . 0579km 

e= . 00093783897 y=2292.2782 

i = 28 . 799652 z = 0. 49842963 

g=358. 27942 Dx=-2 . 2769229 

h=19. 77797 Dy=6. 3292104 

1=1.7260931 Dz=3. 6977685 

ODay 2hr Omin Osec 
a=6774 . 0842 x=-386 1 . 2538km 

e=. 00079610195 y=4620.3061 

i=28. 768322 z=3092.0776 

g=122. 68668 Dx=-6 . 1 86594 1 

h=19. 164909 Dy=-4. 393752 

1=345.68131 Dz=- 1 . 1635444 

ODay 4hr Omin Osec 
a=6776 .5215 x=-4009 . 825km 

e=. 0012429569 y=-5092.6099 

i=28. 787234 z=-1947.7689 

g=215. 66958 Dx=6. 0722793 

h=18. 629512 Dy = - 3 . 6 47 4 1 98 

1=1.0271624 Dz=-2 . 9649566 

ODay 6hr Omin Osec 
a=6776 .6274 x=6330 . 4385km 

e= .00089319126 y=-1514.3192 

i = 28 . 788265 z = -1864.65 

g=307. 46673 Dx=2. 4748519 

h=17. 969212 Dy=6. 6032506 

1=17.622165 Dz=3. 0318842 

ODay 8hr Omin Osec 
a=6773 . 8994 x= 1 4 1 . 965 1 3km 

e=. 00067130355 y=6005.1166 

i=28. 767706 z=3122.1801 

g=97. 726848 Dx=-7. 593806 

h=17. 431316 Dy=-0. 371 13453 

1=335.70208 Dz=l. 0545194 

ODay lOhr Omin Osec 
a = 6778 . 0255 x = -6 42 1 . 5 1 1 5k m 

e=. 0014397351 y=-2136.1552 

i=28. 799631 z=-102.4354 

g=182. 92627 Dx=2. 1778073 

h=16. 822274 Dy = -6 . 3687 1 1 7 

1=358.87728 Dz=- 3 . 6978363 

ODay 12hr Omin Osec 
a=6774 . 02 1 x=3800 . 4452km 

e=. 00073639331 y=-4694.1259 

i =28 . 769 z=-3057 . 2369 

g=267. 80372 Dx=6. 2692341 

h=16. 204838 Dy=4. 2416794 

1=22.382063 Dz=l. 2757673 

ODay 1 4hr Omin Osec 
a=6776 . 1906 x =4 1 03 . 0 1 36km 

e=. 00085551679 y=4988.4883 

i=28. 78614 z=2030 .0187 

g=52. 977531 Dx=-6 . 02 47652 

h=15. 671801 Dy=3. 7762921 

1=345.55649 Dz=2. 8918973 
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Appendix B 




ORIGINAL PAGE IS 
OE POOR QUALITY 


ND 


1 >SAVE "BqDCseculr"LINE 2:PRINT "VERIFY" : VERIFY “ " 

2 DEF FN 1(S,C,T>=(PI AND C<0)+(PI2 AND C>0 AND 3<0)+(ATN T A 
C<>0)+(PI/2 AND C=0 AND S>0)+(3*PI/2 AND C=0 AND S<0) 

3 DEF FN M(X,N)=X-N*INT (X/N) 

4 DEF FN A( A ) =FN M(A,PI2) 

5 DEF FN T(S,C)=(PI2 AND S<0)+(-l AND S<0 OR +1 AND S>-0)*ACS 


C 

6 DEF FN I ( X ) = I NT ( X* 1 E9 + 0 . 5 ) / 1 E9 

8 PRINT DC-ORBIT "’’’'"LOAD ’’OBS ARRAY’’ DATA 0()": STOP 
: DIM 0(6,101): GO SUB 9990 

9 GO SUB 9900 

10 CLS : LET N=6: PRINT AT 20,0;“ solve for DOTa ? PRINT 

AT 20,26;”L“: INPUT D$ : IF D$=“yes” THEN LET N=7 

1 1 CLS : DIM K ( 7 ) 

12 PRINT "INPUT FROM KEYBOARD?": INPUT K$ : PRINT K$ 

13 IF K $ = " n o ” THEN LET al=0(l,l): LET el=0(2,l): LET i I =0 (3,1 

): LET gl =0( 4,1): LET hl=0<5,l): LET 11=0(6,1): GO SUB 1000: G 

O TO 20 

14 IF KS = “ ye s " THEN PRINT "INPUT a I" : INPUT al : PRINT al: PRI 

NT "INPUT el": INPUT el: PRINT el: PRINT "INPUT il”: INPUT il: P 

RINT il: PRINT "INPUT gl": INPUT gl : PRINT gl : PRINT "INPUT hi": 

INPUT hi: PRINT hi: PRINT "INPUT II": INPUT II: PRINT II 

15 LET i I = i I *DTR : LET gl =gl *DTR : LET hi =hl *DTR : LET 1I=1I*DTR 

16 PRINT "INPUT DOTa": INPUT DOTaO : PRINT DOTaO 

20 LET a0=al : LET e0=el: LET i 0= II: LET g0=gl : LET h0=hl : LET 
10=11: LET DOTa=DOTaO 

21 LET k ( 1 ) =aO : LET k(2)=eO: LET k(3) = iO: LET k(4)=gO: LET k(5 
)=hO: LET k(6)=10: LET k(7)=D0Ta 

22 LPRINT •' BgDCsecu 1 r ": LPRINT : LPRINT " INPUT “ 

23 FOR i = 1 TO 6 

24 IF i < = 2 THEN LPRINT " kO( i ;") = ” i k ( i ) 

25 IF i > 2 THEN LPRINT " kO ( “ ; i ; " ) = “ i k ( i ) *RTD 

26 NEXT i 

27 LPRINT "DOTaO=" ; DOTaO 

28 PRINT " GOTO 8000 to Test BgDCseculr " 

29 STOP 

30 REM DIFFERENTIAL-CORRECTION 

31 CLS 

32 LET C$="NO" 

34 LET I T=0 

37 GO SUB 500 

38 DIM M ( 7 ) : DIM D(7, 101 ) 

39 REM A . dK=dO 

40 FOR t =0 TO ENDt STEP DELt 

41 LET nt=t/DELt+l 

42 CLS : PRINT nt 

44 GO SUB 270 

45 GO SUB 400 

50 FOR 1=1 TO N 

51 LET B(I)=D(I,nt)+B(I) 

52 FOR J=1 TO N 

53 LET A( I , J)=S( I , J)+A( I , J) 

54 NEXT J: NEXT I: NEXT t 

55 REM NORMAL EQ . LINEAR IN Dn 

56 GO SUB 505: CLS : BEEP .5,0: BEEP 1,12 

57 FOR i = 1 TO N: LET k ( i ) = k ( i ) -t-C ( i ) : NEXT i 

58 LET a0=k ( 1 ) : LET eO=k(2): LET iO=k(3): LET g0=k(4>: LET hO= 
k ( 5 ) : LET 1 0 = k ( 6 ) : IF N = 7 THEN LET D0Ta = k(7) 

59 PRINT " parameters DP " ; TAB 25;"# ";IT: LPRINT ” mean e 

lements“;TAB 20;'" DEL.ele.# ";IT 

60 FOR 1=1 TO 6 

6 1 IF I < = 2 THEN PRINT ” K ("; I ;”) = "; FN I(K(I)),FN I(C(I)) 


62 IF I < =2 THEN LPRINT M K ( H ; I ; " ) = " ; FN I(K<I));TAB 20;FN I(C(I 

) ) 

63 IF I > 2 THEN PRINT " K < " ; I ; “ ) = " ; FN I ( K ( I ) ) *RTD , FN I(C(I))*RT 
D 

64 IF I > 2 THEN LPRINT " K ( " ; I ; “ ) = “ ; FN I ( K ( I ) ) *RTD ; TAB 20;FN I( 
C ( I ) ) *RTD 

65 NEXT I: IF N = 7 THEN PRINT "K<7)=";FN I(K(7)),FN I(C(7)>: I 
F N=7 THEN LPRINT " K ( 7 ) = M ; K C 7 ) ; TAB 20;C<7) 

66 IF ABS C( 1 ) > 1 E-5 THEN GO TO 68 

67 LET C$ = " YES" : GO TO 70 

68 LET I T= I T+ 1 

69 IF I T> = 1 1 THEN GO TO 88 

70 REM AVE. REL . ERROR 

71 FOR K=1 TO 6 

72 LET SUME=0 

73 FOR 1=1 TO MAX 

74 IF OC K , I ) < >0 THEN LET RELE = ABS D( K, I ) /0( K , I ) 

75 LET SUME=SUME+RELE : NEXT I 

76 LET AVE=SUME/MAx 

77 PRINT " ( " ;k; " )-ERR= * ; AVE*100; "%" 

78 LPRINT * ( M ;K; “ )-ERR= M ; AVE*100; 

79 NEXT K 

80 LPRINT : CLS : GO TO 87 

81 FOR i=t TO n 

82 PRINT i 

83 FOR k=l TO max 

84 IF i < = 2 THEN PRINT k,D(i,k> 

85 IF i > 2 THEN PRINT k,D(i,k)*RTD 

86 NEXT k: NEXT i 

87 IF C$ = " NO" THEN GO TO 37 

89 PRINT : PRINT " conv.ele. M 

90 LPRINT : LPRINT " conv.eie. " 

91 FOR k=l TO 6 

92 IF k< = 2 THEN PRINT " k < * J k ; * ) = " ; k ( k ) 

93 IF k< = 2 THEN LPRINT H k ( " ; k k ( k ) 

94 IF k > 2 THEN PRINT " k ( " ; k ; " ) = " ; k < k ) *RTD 

95 IF k> 2 THEN LPRINT " k ( " ; k ; " ) = " ; k ( k ) *RTD 

96 NEXT k 

97 PRINT "DOTa=“ ; DOTa 

98 LPRINT DOTa= " ; DOTa 

99 STOP 

200 REM BgSecul arORBI T 

201 LET ADP=aO : LET EDP=e0: LET IDP=i0: LET GDP=gO: LET HDP=hO: 
LET LDP= 1 0 

202 LET NO=SQR (GM/ADP~N3) 

203 LET EDP2=EDP*EDP : LET CN2=N1-EDP2: LET CN=SQR (CN2 ) 

204 LET GM2=K2 / ADP / "N2 : LET GMP2=GM2 / ( CN2*CN2 ) : LET GM4=K4/ADP*N 
4: LET GMP4=GM4/CN^8: LET F 1 D4G2 = F 1 D4*GMP2 

205 IF Dt =0 THEN LET CI=COS (IDP): LET CI2=CI*CI: LET CI3=CI2* 
Cl : LET C I 4 = C I 2*C I 2 

206 REM 1 DOT , gDOT , hDOT 

207 LET 1 DOT = NO* ( CN* ( GMP2* ( F3D2* ( N3*CI2-N1) +F3D32 *GMP2*( N25*CN2 
+ 1 6*CN- l 5 + ( N30-96 *CN-N90*CN2 )*CI2+(N25*CN2+1 44*CN+ 1 05 ) *CI 4) ) )+Fl 
5D16*GMP4*CN*EDP2*<N3-N30*CI2+N35*CI 4) ) 

208 LET nD0T=N0+ 1 DOT : LET PD = PI2/nD0T 

209 LET gD0T=N0*< (GMP2* ( F3D2*( N5*CI 2-N 1 > +F3D32*GMP2*< N25*CN2+24 
*CN-N35+(N90- 192*CN- 126*CN2)*CI 2+< 45*CN2+360*CN+N385 >*CI 4) ) ) + F5D 
1 6*GMP4*( 2l-N9*CN2+(N126*CN2-270)*CI2+(N385“189*CN2)*CI 4) ) 

2 10 LET hDOT = NO* < ( GMP 2* < F3D8*GMP2* ( (N9*CN2+12*CN-N5)*CI-<N5*CN2 
+36*CN+N35)*CI 3>-N3*CI ) ) +F5D4*GMP4* ( N5-N3*CN2 ) *CI *( N3-7*CI 2 ) ) 

211 RETURN 

220 REM SECULAR -GDP , HDP , LDP 
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S G G 0 D P p ;?n<CD*PK IP =OP<0 THEN LET COP-GEP^.2 
tf T S&rXZSE IP hdp<0 then let hdp.hdp*pi 2 

225 LET L0P»1°*P“°I*“! ,p LDP <0 THEN LET h-HDP s lET 

*“ g? SS LET eUop! LET LET .-«*• 

1 =LDP 

229 RETURN 

270 REM ORBGEN x H EN GO SUB 200 

272 IF T D0Ta=0 THEN GO TO 290 Q; G0 T0 289 

| 73 XF t=0 THEN LET LE ^ THEN let slgD .-l 

! sr» s " ,ot,pD> 

1 1 Spfi|E-T S,SS3f!STSa h 

282 LET aO=aO+DELa : LET eu e 
^B^LE^DfPDs LET tsum = tsum+PD 

HI GO SUB 200: GO SUB 220 LET HO=H: LET 10-1 

285 LET aO=a : LET eu e. 

287 NEXT k 

288 GO SUB 200 

289 LET tO=tsum : LET Dt- m(5 )=h: 

250 GO SUB 220 «C3»-1= LET " l4> - 9: LET 

295 LET m(l)=a: LET mC2J 

LET m(6) = l , n / to) 

299 RETURN State matrizart O 

300 REM Secular retarded 

Si Li” B2=C2,p2 

306 LET F^A=l/a 0 : LET FEDP=e0/p LET S 3 = 2 *B 2 *n 0 *SIN 10: LE 

307 LET S 1 =F7D2*F 1 DA : LET SC 

T S4=F3D2*F lDA*nO 

308 LET SO, 0 = 1 

309 LET SC 2 , 2 > = 1 

310 LET S( 3, 3) = 1 

31 1 LET SC 4 , 4> = 1 

312 LET SC 5 , 5 ) = 1 

313 LET SC6»6)-1 

319 RETURN 

320 REM t 

III So!n = l -S^*9mD0T*t 

“ S *<*.*>- 

324 LET SC 4 , 3 ) = -5*S3*CI*t 

S K? s<5!n=i a it.h.D0T,t 

327 LET S<5,2>- 4*S2*h«D0T»t 

328 LET SC 5 , 3 ) = S3*t 

329 REM 1 part ial s nT v* t 

^ r ut eft i > = -CS4 + Sl *lmDOT 

332 LET S<6.3>- - 3«cn«S3.CI 



340 RETURN 


350 REM 

352 LET 
A=n0/a0 

353 LET 
03*F 1 DA 

355 LET 

369 

370 REM 

371 LET 

372 LET 

373 LET 


DOTa 

F1DA= 1 /aO: 
D4=( 1-eO) ; 


LET pO-aO*( 1 -eO*eO) 
LET D3=FEDP*D4 


374 

375 

378 

379 

380 

382 

383 

384 

386 

387 

388 

399 

400 
402 

> 0 ) 

405 

406 

408 

409 

410 
412 


S( 7 , 7 ) = 1 
RETURN 

retarded 
t2=t*t 
S( 1 ,7> = t 
S(2, 1 ) = 

S ( 2 , 2 ) = 

S ( 2 , 7 ) = 


: LET FEDP=eO/pO: 
LET D2=DOTa*eO: LET 


part ial s 


LET 

LET 

LET 

LET 

LET 

LET 

LET 

LET 

LET 

LET 

LET 


04*F l DA*F 1 DA*DOTa*t 
1 -DOTa *F 1 DA* t 
D4*F 1 DA*t 

1 > ~ 4 *Dl*gmDOT*t2 
1 4 ’? 4 ' 2 )+ 2 * D 2 *9">D0T*t2 

| ( 5 # [ Jiff 4 ’ J>*t+4*D3*gmDOT*t2 
1 )_4 *Dl*hmDOT*t2 

q^'7! = f <5 ' 2)+2 * D2 * hmD0T *t2 

S(6 ' 1) -ttl’ ! >*t+4*D3*hmDOT*t2 

s < 6 : 2 ) : s < 6 ,' 2 ): o 2 » i 4 ; Do ? 2 * : '' ,DS ‘ D " M * DOTa * t2 


TFT o / £ n \ rs , - - - 1 miuu i 

RETURN ” 6 ’ 1 ) * t-F 3D2*FlD2*DNDA*t2+3*D3*lmDOT*t2 


REM part ial s 
DEF FN D(X, Y)=( PI 2-ABS 


FOR 
LET 

IF k<=2 THEN 
IF k>2 THEN 
NEXT k 

IF ABS D(4,nt)>Pl THEN 
(4 ;" t ^ FN D <x,y) : PRINT "D< 
414 IF ABS D( 6 , nt ) > PI THEN 


X >*< 1 AND Y<0)+(Pi 2-ABS 


k= 1 TO 6 

D(k,nt)= 0 (k,nt) - m ( k > 

PRINT "k=";k,D(k,nt) 
RINT "k=" ;k,DCk,nt)*RTD 


X>*< 


( 6 , at ) -FN D( x , y ) 
425 IF t=0 THEN 
GO SUB 320 
IF N=6 THEN 
IF t=0 THEN 
LET D0=(0(1, 
LET DC= ( n> ( 1 ) 


426 

430 

432 

436 

438 

440 

445 

470 

490 

495 

499 

500 
502 

504 

505 

506 

507 

508 

509 
51 1 
512 
5 1 4 
515 
517 
5 18 


PRINT 
GO SUB 


D< ” 
300 


. kET x=D(4,nt): LET 
* 4; ) = * ; D< 4, nt ) *RTD 
LET x=D(6,nt): LET 
'6;")=";D(6,nt) *RTD 


y=D(6,nt: 
y=D< 4 , nt 1 


return 

go TO 470 
nt ) -0< 1 ,nt- 
mO) /DELt 


1 > ) /DELt 


350 


LET D(7,nt)=D0-DC 
PRINT ”k=7",D(7,nt) 

LET mO=m ( i ) 

IF t=0 THEN GO SUB 
GO SUB 370 
RETURN 

REM GAUSS elimination 

DIM A < 7 , 7 ) : DIM B(7)- 
RETURN 

REM PIVOTAL CONDENSATION 
GO SUB 5000 
NM 1 =N- 1 
K=1 TO NM 1 
kPl=K+l : LET L=K 
I=KP1 TO N 

A( I , K ) >ABS ACL, K) 


DIM C ( 7 ) 


REM 
LET 
FOR 
LET 
FOR 
IF ABS 
NEXT I 
IF L = K THEM 
FOR J=k TO N 
LET S=A(K, J) 


THEN LET L=I 


GO TO 524 


LET DND 
D1 =DOTa* 


AND Y 


LET D 
LET D 
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NEXT J 
NEXT K 


519 LET A(K, J>=A(L, J> 

520 LET A(L,J)=S 

521 NEXT J _ . r ft rh )=S 

522 LET S=B(K) : LET BCK)=BCL): LE 

523 REM ELIMINATION 

524 FOR I=KP1 TO N 

525 LET PI VMPY= A< I , K) / A ( K , K ) 

526 LET A( I ,K)=0 
K 07 FOR J=KPl TO N 

528 LET A<I,J>=A(I,J> - J J 

529 LET B(I)=B(I)-PIVMPY*B<K) : NEXT I 

530 LET C(N)=B(N) / At N , N> 

1 E »hJs»” next j 

536 LET CC I )=(B( I )-SUM) /A< I , I > : NEXT L 
550 RETURN 

1000 REM calculate DOTaO 

1005 LET sum=0 _ . . . . n , . ._i). NEXT i 

1010 FOR j = 2 TO max: LET sum = sum + 0< 1 , j 1 0(1, j 1 >- ne. 3 

1020 LET DQTaO=sum/ENDt 

1099 RETURN 

5000 REM MATRIZART 

5001 CLS 

5005 PRINT “MATRIZART" 

5006 LPRINT "MATRIZART" 

5010 FOR 1=1 TO N 

5020 PRINT ?B(I> 

5021 LPRINT "B(";l;“)="iB(D 
5030 FOR J=1 TO N 

* ’ ) = " ; A< I , J > : 

“ ) = " ; A ( I , J ) : 


A ( 
A( 


I ; J; ‘ 
1 } I i J J 


5040 PRINT 

5041 LPRINT 
5050 NEXT J 

5060 PRINT 

5061 LPRINT 

5065 PAUSE 300 

5066 CLS 
5070 NEXT I 
5090 RETURN 

8000 REM Test BgDCmean 

8001 INPUT " DELt : " 5 DELt : INPUT 

8002 LET MAX= ( ENDt /DELt ) + 1 
8005 FOR t=0 TO ENDt STEP DELt 
8010 LET nt=t/DELt+l 
8020 GO SUB 270 
8025 CLS 

8030 LET 0( 1 , n t ) = a : LET 0(2,nt)=e 

8035 FOR* K= 1 = T0 6 V NT* C i<K . nt > AND K< = 2 , + OCK , nt > *(RTD AND K>2 

) ’ : NEXT K 
8040 NEXT t 

8050 BEEP .5,0: BEEP 1, 12 

8099 STOP : GO TO 9 

9900 REM "CONSTANTS" 

9901 RESTORE 

9920 DATA 398600^63! 6 378^ i 66 ! 0^7292 1159E-4, 298. 25,-0. 1 08265 1 7E- 2 

9924 DATA . 5 , . 3 3 3333 3 333 , . 25 , 1 . 5 , . 37 5 , . 09 37 5 , 1 . 25 , . 3 1 25 , .9375,3. 
9926 READ N 1 , N2 , N3 , N4 , N5 , N6 , N9 , N 1 1 , N25 , N30 , N35 , N40 , N90 , N 1 26 , N385 


• ENDt : * 5 ENDt 


LET 0<3,nt)=i: LET 0(4,nt)-g: 
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^oL^Ki;^-^;^^! 1 - 25,30 - 35 ' 40 - 90 - 126 ’ 385 

9940 LET K2=-0.5*J2*Re~2 
9945 LET C2=F3D2*K2 
9950 LET K3=J3*Re~3 
9960 LET K4 = F3D8*J4*Re , '4 
9970 LET PI 2=2*PI 

9980 LET RTD= 180/PI : LET DTR=PI/180 

9989 RETURN 

9990 REM OBS. ARRAY 

9991 LOAD "OBS ARRAY" 

9992 FOR 1=1 TO 100 

9993 PRINT 

9994 FOR K= 1 TO 6: 

: NEXT K 

9995 NEXT I 

9997 PRINT E "ni?f^ LET ENDt=0< 2 » 1 0 1 ) : LET MAX= (ENDt /DELt > + 1 

9998 RETURN > DELt ’ "ENDt ="; ENDt '" MAX= MAX 

9999 STOP 


DATA 0( ) 


PRINT (OIK, I) AND K< = 2 ) +0 ( K , I ) *(RTD AND K>2)' 
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1& BgDCseculr 


INPUT 

kOCl>=6775.8712 
kOC 2 ) = . 0005736 1504 
k0( 3)=28 .782549 
k0( 4 ) = 247 .80866 
k0(5)= 19 .779807 
k0 C 6 ) = 1 12.19 167 
DOTaO=-4 . 6642327E-6 
BgDCsecu I r 


INPUT 

kOC 1 >=6778 . 1 4 
kOC 2)= . 001 
k0(3)=28 .8 
k0( 4) = 360 
kO(5)=19.78 
k0<6)=0 
DOTaO=0 
mean elements 
KC 1 ) =6775 . 8586 
K(2>=. 00057168 
KC3>=28 .782549 
K<4)=314. 49322 
K<5)=19. 781771 
KC6)=40. 79636 
K<7)=-4. 4657805E-6 

< 1)-ERR=. 042215133% 
(2)-ERR= 108 .07809% 

< 3>-ERR= . 060630488% 

(4)-ERR=42. 376957% 

( 5 >-ERR=0 . 7 1 432865% 
( 6 ) -ERR = 1 13.89555% 


DEL.ele.# 0 
-2 . 2814093 
- . 00042832 
-.017450977 
-45 . 506783 
.0017712417 
40.79636 
-4. 4657805E-6 


mean elements 
K( 1 ) =6775 . 8707 
K<2>=. 000573535 
KC 3) =28 .782549 
K(4>=247. 80838 
K<5)=19. 779979 
KC 6 > = 1 12.1 3376 
K(7)=-4.6557707E-6 
( 1 >-ERR= . 0002384396 
<2)-ERR=0. 52303362% 
( 3)-ERR=0% 

( 4 ) -ERR=25 . 195103% 

( 5 ) -ERR= .051707417% 

(6) -ERR=68 . 978533% 


DEL.ele.# 1 
.012070939 
1 . 855E-6 
0 

-66.684835 
-.0017924412 
7 1 . 337405 
- 1 . 8999025E-7 
3% 


mean elements 
K ( 1 ) =6775 . 87 1 2 
K(2)=. 00057361 1 
KC 3> = 28 .782549 
K( 4)=247 . 80852 
K(5)= 19 .779896 
K(6)= 1 12. 19086 
K(7)=-4.6638612E-6 


DEL.ele.# 2 
.000514262 
7 . 7E-8 
0 

. 00013464508 
-.000083136176 
.057097192 
-8 . 0904803E- 9 


( 1 ) -ERR= . 0000 10145558% 

(2) -ERR=. 02255187% 

( 3 ) -ERR=0% 

( 4) -ERR= .000055049224% 
C5>-ERR=. 001 202 1346% 
(6)-ERR=. 058431 165% 


mean elements 
KC 1 >=6775.8712 
KC2)=. 000573615 
KC3)=28. 782549 
KC4)=247. 80852 
KC5>=19. 779892 
K C 6 ) = 1 12. 19328 
KC7)=-4.6642l2E-6 
C 1 ) -ERR = 4 .6109305E-7% 

C 2 > -ERR= . 000977 14713% 

C 3)-ERR = 0% 

C 4)-ERR= .000026798669% 
C 5 > -ERR= . 0014324659% 
C6)-ERR= .0026780289% 


DEL.ele.# 3 
.000019898 
3E-9 
0 

5 . 6722822E-6 
-3 . 5523383E-6 
.0024229239 
-3 . 508 2625E- 1 0 


mean elements 
KC 1 ) =6775 . 87 1 2 
KC2)= .000573615 
KC 3)=28 .782549 
KC4>=247. 80852 
KC5>=19. 779892 
KC 6 ) = 1 12. 19339 
KC7)=-4.6642192E-6 
C 1 > -ERR = 2 . 1 303428E-8% 
C2)-ERR=. 000048678 18 2% 
C 3>-ERR=0% 

C 4) -ERR=. 000026753357% 
C 5 > -ERR= . 00 1 506 1 207% 

( fi)-FRR= .00092556834% 


DEL.ele.# 4 
2 . 37 IE-6 
0 


0 

2 . 864789E-7 
- 1 . 7 1 88734E-7 
.00010777336 
-7 . 1597023E- 1 2 


conv .ele . 
kC 1 >=6775.8712 
kC2>=. 00057361481 
kC3)=28. 782549 
kC4)=247. 80852 
kC5)=19. 779892 
k C 6 ) = 1 12. 19339 
DOTa=-4 . 66421 9 2E-6 




Appendix C 



1 >SAVE -BgOBSWEAN-LINE 2.PRI»T ;VEEirV- : »ERIFY D -- (o)t(JTN T „ 

ND l^^ S ^ T ^0 P ^>0 *» c - 3<0> 

S 

5 DEE FN R<X.Y.Z>.SOR <X.X*Y*YtZ*Z> 

I E H AND S<0 OF »ND S>- 0 ).ACS 

’ .! |;rSSssi:5.r.;siK's;»: 

u PrInT -INPUT- K -FOR Kepler ian elements inpu 
1 3 4 PRlS MNPUT FROM KEYBOARD INPUT KS 

il “"^aSS&s a=T;; m , 

18 IF 0 $= M C“ THEN GO TO 29 

II ?^»T N -I T Np5f“:‘?»P^ A, PRINT a 

21 PRINT “INPUT e"r INPUT e: P®™* ^ i=i * D TR 

22 PRINT “INPUT 1 : INPUT 1 . P _ , *htd 


PRINT e 
PRINT i 


LET i=i*DTR 


23 PRINT 

24 PRINT 

25 PRINT 


INPUT g" 
INPUT h" 
INPUT 1” 


InPUT g: PRINT g: LET g=g*DTR 
INPUT h: PRINT h: LET ^=h*DTR 
INPUT I: PRINT 1 ’■ LET 1 I D 


z “ r **" 1 a...-- PR T N T 1: LEI l=i*uir 

25 PRINT -INPUT I s INPUT 1. P R™T , . /DTR .- h = 

26 LPRINT “ a = 5 a e = -e 1- ' 1/D1K 9 

1 / D TR 

27 GO SUB 300: GO TO 40 

28 REM INPUT CARTESIAN »»CUL’’“: INPUT I 

29 PRINT "INPUT UNITS ''KM 0* <=UL 

30 PRINT "INPUT x : INPUT *. 


; h/DTR' ” 1 = ' 


INPUT U$ 


31 PRINT “INPUT y“: 

32 PRINT "INPUT z" : 

33 PRINT "INPUT Dx" 

34 PRINT "INPUT Dy _ 

35 PRINT “INPUT Dz“ 

36 PRINT " x = ” J x ’ " y = 


INPUT y: PRINT y 
INPUT z: PRINT z 
INPUT Dx: PRINT Dx 
INPUT Dy: PRINT Dy 
INPUT Dz: PRINT Dz 


INPUT Dz“: I^ Ul .. U :;.nx- : ;Dx'-Dy = -;Dy--Dz=-;Dz 
x = ” ; x ’ " y = ;y : z _, .. nv= » ; n *’ »dv=“ ;D y* “Dz =“ ;dz 


36 PR1N1 X= 7- -nx=" :Dx’ “Dy= »Dy uz= 

ll xr R ;"--cJl : ^HEN = let X=X*Re : LET y = y*Re: LET z = z*Re : LET D 
x=Dx*Re*Ke : LET Dy = Dy*Re*Ke: LET Dz=Dz*Re *Ke 

39 GO SUB 400 

« *“ LET K<3) = i : LET KM>-» = LET K<5.=h: 

Tz LET ^ X ( 1 > = X : LET X<2>-y: LET X.3« LET X<4..Dx: LET X< 5 >-D 

y: LET X(6)=Dz 

45 REM CALL OSCMEAN 

46 GO SUB 1000 

47 GO SUB 4000 

48 STOP : CLEAR : GO TO 2 

50 REM INPUT DATA ARRAY: GO SUB 9995 

51 FOR t = 0 TO ENDt STEP DELt 

52 LET nt=t/DELt+l y =P(l,nt,2): LET z=P<l 

„;y,T 00 

THEN LET JHP.2.J. •>>; «t e - E P<Z, ? t.2>^LET ,=P<2 
,nt , 3) : LET g = P(2,nt,4): LET h-P<2,nt . 0 J . 

B 300 


C-l 


55 

LET 


LET K < l)=a: 
K < 6 ) = 1 : LET 


LET K( 2 ) =e : 
X(l)=x: LET 


X( 5 ) =Dy : LET X(6)=Dz 
56 REM CALL OSMEAN 
GO SUB 80 
GO SUB 1000 
LET G<nt , 1 ) =ADP 
, 4)=GDP*RTD: LET 
GO SUB 4000 
NEXT t 

BEEP .5,0: BEEP 
STOP : CLS 
REM DOTa 
LET D0Ta=0 

FOR t=0 TO ENDt STEP DELt 
LET n t = t /DELt + 1 

IF nt > 1 THEN LET dela=G(nt 1 ) 
LET DOTa=DOTa+de la ' 

NEXT t 

LET nDEL=(nt- 1 )*DELt 
LET DOTa=DOTa/nDEL 
PRINT BRIGHT i ; AT 11 


LET K( 3) = i ; 
X ( 2 ) =y : LET 


LET K ( 4 ) =g ; 
X(3)=z: LET 


LET K(5) 
X( 4) =Dx : 


: h : 
LE 


57 

58 

59 
G < n t 

60 
66 

67 

68 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 
90 

80 
81 
82 

83 

84 
c+.5) 

8 6 CLS 
87 

89 

90 

92 

93 

94 

95 

96 


‘ EDP; LET G ( nt , 3 ) = I DP*RTD • LFT 

G<nt,5)=HDP*RTD: LET G C n t , 6 ) =LDP*RTD 

1,12 


-G(nt- 1 , l ) 


REM 

LET 

LET 

LET 

LET 


,4, DOTa- DOTa km/sec " : STOP 

t — Da y , hr , min, sec 
Day= I NT (t/86400) 
hr = I NT ( ( t-Day*86400> /3600) 
m l n= I NT ( ( t-<Day*86400+hr*3600> ) /60) 
sec=INT Ct-(Day*86400+hr*3600+min*60)) : 


GO TO 


INVERSE 1 ; Day ; "Day ";hr;"hr 


elements 
STEP DELt 


PRINT 
RETURN 

REM output mean 
FOR t=0 TO ENDt 
LET nt=t/DELt+l 
GO SUB 80 

PRINT "MEAN ELEMENTS 
PRINT -a"=“;G(nt,l) 

< nt , 4 ) ' “h ' ’ = “;G(nt,5) ' "l ' 

97 PAUSE 4E4 

98 NEXT t 

99 STOP 

100 REM SOLVE KEPLERS EQ 
1 10 LET EA=0 

115 IF 1=0 THEN GO TO 160 
120 LET EA= 1 +e 

125 FOR N=i T0 10 . LET 0EA=EA: L£t EE=EA _ e * SIN 
FE/ ( 1 -e *COS ( EA-0 . 5*FE ) ) : LET DEA = ABS (EA-OEA) 


LET sec=INT <se 


imin;“min ”;sec;"sec‘ 


'“e” = ”;G(nt, 2 ) 
= J G ( n t , 6 ) 


i ' ’ = " ; G(nt , 3) 


135 
1 40 
160 

199 

200 
201 
LET 
202 

203 

204 


EA-1: LET EA=EA- 


■ IE-8 THEN GO TO 160 


EDP-eO: LET IDP=iO: LET 


GDP=gO : LET HDP=hO: 


IF DEA< =0 . 

NEXT N 

LET EA=FN M(EA,2*PI) 

return 

Rem BgORBIT 
LET ADP=aO: LET 
LDP = 1 0 

LET N0=SQR (GM/ADP / 'N3) 

“20^^! 1 = S ^ T ' 4 " 

cf 05 Jr F SiJcS:® CI ’ C0S «T CJ2.CI.CI: 

230 REM SP , LP- CONSTANTS 


LET Cl 3=CI 2* 
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LET S2G=N2*SG*CG: 


LET H 


in 2 m L wrSf'r C !fr N: LET CN6 = CN3 *CN3: LET F 1D1 CN= 1/(1 +CN) : LET F 
1DCN3=1/CN3: LET F1DCN6=1/CN6 

233 LET GM3 = K3/ADP' V 3: LET GMP3=GM3/CN6 : LET G3DG2=GMP3/GMP2 

234 IF Dt-0 THEN LET SI=SIN (IDP): LET TI=SI/CI: LET P3T2M1=N3 

LET P 1 MT2 = N 1 -Cl 2 : LET SQ1MT2=SQR (P1MT2): LET T31MT2=N3 
* P ^ 2 = LET T5T2M 1 =N5*CI 2-N 1 : LET P3M5T2=N3-N5*CI 2 : LET A0=CI2/( 
N1 N5*CI 2 ) : LET A1=FID2*F1D4*(N1-N1 1*CI2-N40*CI 2*A0) • LET A3--F1 
D2*F1D4*CI*(N1 1 +80* A0+ 2 00* AO* AO) LET A3- FI 

P6=?i*SP3*SqIm^ 2 3 * EDP: LET SP3=F1D2 * GMP2: Lp T TSP3=CI*SP3: LET S 

236 LET A2=CN3*GMP2*Al-FlD4*FlD4G2*(N2+EDP2-400*EDP2*CI2*AO*An- 

240 REM UDP, PERIODIC TERMS 

P^FN A E UDP P=EDP: LET GP=GDP: LET LP=LDP: LET UDP=GDP+LDP: LET UD 
242 REM LP-TERMS 

2G=N2*CG*CG-N 1 ^ (GDP>: LET CG=C ° S <GDP>: LET S2 G=N2*SG*CG : LET C 

ip 4 ^».o ET D 1 E=A4*SG+EDP*GMP2*A 1 *C2G : LET D 1 1 =- (EDP*D IE) / TI • TFT n 
D2E ' EDP ‘“3*GM P 2,AUS 2 G-CN3.A4,M L " ° 

245 LET EP=SQR (D2E*D2E+ (EDP+D 1 E) * C EDP+D 1 E) ) 

P=FN A(Jp) P=HDP+EDP2 * A3 * GMP2 * S2G+(( EDP*CI*A4)/(SI*SI))*CG: LET H 
upIfn E I(up) UDP+A 2 * S2G+<(EDP * A4 * F1d1CN) * (N2+cn_EDP2)+A5) * CG: let 

248 LET SL=SI N (LDP): LET CL=COS (LDP) 

249 IF EDP> =0.05 THEN LET SM=D2E*CL+(EDP+DtE)*SL- LET CM-(Fn P + 

(GP): LET S2G=N2*SG*CG : LET C2G=CG*CG-N1 C ° S 

250 REM FP 

)f 5 LET" E L"c05 <eL*' KP! 00 SUB 100: “ T EAP=ES: LET SEA-SIN <EA 

253 REM SP-TERMS 

254 LET CF2 = CF*CF : LET CF3=CF2*CF: LET S2F=N2*SF*CF • LFT r?F-M-> 

^dp 2 o N1: LET S3F=N3 * sf_N4 *SF*SF*SF: LET C3F = N4*CF3-N3*CF • LET 
CPF.S2G.CF.C2G.SF: LET S2GP2F.S2G.C2F+C2G.S2F* LET 

LET 

255 REM COMPUTE a,e,i,g,h,l 

^256 LET a=ADP* (N 1 +GM 2 * ( P3T2M 1 *( ADR3-F 1DCN3)+T31 MT2*ADR3*C2GP2F) 

s:|ipirrsiiSSisSsS 

=FN 8 A(n 1=IDP+D1I+SP6 *<N3*C2GP2F+EDPT3*C2GPF+EDP*C2GP3F): LET i 
*S2GP3F )" ) ^ = LET T h = FN ^ A ( h) FP_LP+EDP *^ F * ~ ^3*S2GP2F+EDPT3*S2GPF+EDP 

CM 2 <>^ L ?H T EN M= LE 2 r?M:SM%M +D1E) * SL: ^ CM =' < EEP ^ E) *CL-D2E*SL : IF 


LET SEA=SI N 
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EA 


g 

h 


262 LET 1=FN 1 <SM, CM, TM) 

264 LET g=u-l: LET g = FN A(g) 

269 RETURN 
300 REM KEP-POSVEL 
302 GO SUB 100 

304 LET cosEA=COS EA : LET sinEA=SIN 
308 LET e l =a*SQR <l-e*e> 

310 LET r=a*( l-e*cosEA) 

312 LET DotEA=SQR (GM/a)/r 
314 LET Xw=a*(cosEA-e ) 

316 LET Yw=e 1 *s inEA 
318 LET DXw=-a*DotEA*s inEA 
320 LET DYw=el*DotEA*cosEA 
322 LET sini=SIN i: LET cosi=COS l 
324 LET s ing=SIN g: LET cosg=COS 

326 LET s i nh=SI N h: LET cosh-COS 

328 LET Px=cosg*cosh-s ing*s inh*cos i 

330 LET Py =cosg*s i nh+s i ng*cosh*cos i 
332 LET Pz = s i ng*s i n i 

334 LET Qx = -sing*cosh-cosg*s inhAcosi 
336 LET Qy=-s i ngAs i nh+cosgAcoshAcos i 

338 LET Qz=cosg*sini 

339 REM x,y,z 

340 LET x=XwaPx+YwAQx 
342 LET y=Xw*Py+Yw*Qy 

344 LET z=Xw*Pz+Yw*Qz 

345 REM Dx , Dy , Dz 

346 LET Dx=DXwaPx+DYwaQx 
348 LET Dy=DXw*Py+DYw*Qy 
350 LET Dz=DXwaPz+DYwaQz 

399 RETURN 

400 REM POSVEL-KEP 

404 LET r=SQR ( x*x+y*y az*z ) 

406 LET V 2 =Dx*Dx+Dy*Dy +Dz*Dz 
408 LET Uz=z/r 

410 LET Wx=y*Dz-z*Dy : LET Wy=z aDx-x aDz . LET 

412 LET W=SQR ( Wx*Wx+Wy*Wy+Wz*Wz ) 

414 LET Wx=Wx/W: LET Wy=Wy/W: LET Wz-Wz/W 

416 LET rDotV=x*Dx+y*Dy+z*Dz 

417 LET Dr=rDotV / r 

418 LET Vz=(r*Dz-Dr*z) /W 

420 REM a 

422 LET a=(GM*r)/(2*GM-r*V2> 

430 REM e 

432 LET Se=rDotV/SQR <GM*a): 

434 LET e =SQR (Se *Se +Ce *Ce ) 


LET Ce= 1-r/a 


440 REM i „ , 

442 LET sin i =SQR < Wx*Wx+Wy*Wy ) 

444 LET cosi=Wz 

446 LET i =FN T(sini.cosi) 

450 REM h ^ . . , 

452 LET sinh^Wx/sini: LET cosh- Wy/sin 

454 LET h=FN T(sinh,cosh) 

460 REM EA _ , 

462 LET sinEA=Se/e: LET cosEA-Ce/e 

464 LET EA = FN T ( s i nEA , cosEA ) 

470 REM f 

47 1 LET EE= 1-Ce 

472 LET s inf =SQR ( 1 -e *e ) *s l nEA /EE 
474 LET cost = CcosEA-e ) /EE 

476 LET f =FN T(sinf.cosf) 

480 REM u 


Wz=x*Dy-y*Dx 
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482 

484 

485 
488 
490 
492 
499 

1000 

1001 

1005 
K( 5 ) 

1006 


LET s i nu=Uz /s i n i : LET cosu=Vz / s i n i : LET cosu=FN I(cosu) 
LET u = FN TC 5 i nu , cosu ) 

REM g 

LET g = u- f : LET g = FN ACg) 

REM 1 

LET 1 =EA-Se : LET 1=FN A < 1 ) 


1 =EA-Se : 

RETURN 
REM OSCMEAN 
LET Dt =0 
LET aO=KC 1 ) : 

LET 1 0 = K ( 6 ) 

LET xM = X C 1 > : 

M = X C 5 ) : LET DzM = X(6) 
1009 LET OLDEL=999 
FOR m= 1 TO 10 
REM CALL BgORBIT 
GO SUB 200 
REM KEP- > POSVEL 
GO SUB 300 
LET ao=a: LET eo=e 


LET e 0 = K ( 2 ) : LET iO=KC3> 
LET y M = X C 2 > : LET zM = X(3) 


LET gO=K( 4) : LET hO= 
LET DxM=XC 4 ) : LET Dy 


1010 

1019 

1020 
1022 

1025 

1026 

1030 

1031 

1032 

1033 

1034 

1035 
1040 
1042 
z = DzM 


LET D< 1 >=X( 1 )-x : 

LET DC2)=XC2)-y: 

LET D(3)=X(3)-z: 

LET D( 4)=X( 4)-Dx 
LET DC 5 > = X C 5 ) -Dy 
LET D(6)=X(6)-Dz 
REM POSVEL- > KEP 
LET x=xM : LET y=yM: 


LET i o= i : LET go = g: LET ho = h: LET lo = l 
D C 1 > 

DC 2) 

DC 3 ) 

+ DC 4) 


LET xM = xM + 
LET yM= yM + 
LET zM=zM + 
LET DxM=DxM 


LET 

LET 


DyM=DyM 

DzM=DzM 


DC 5 ) 
DC 6 ) 


LET z=zM : LET Dx=DxM: LET Dy=DyM : LET D 


1045 GO SUB 400 

1046 LET aO=a : LET eO=e: LET iO=i: LET gO=g: LET h0=h: LET 10=1 

1049 CLS 

1050 PRINT TAB 9; INVERSE 1 INTER .#" J m 

1051 LPRINT : LPRINT TAB 8;”Iteration # : " J m 
1053 IF m> 1 THEN GO TO 1065 

1055 PRINT ” Kepler input Cartesian” 

1056 LPRINT ” Kepler input Cartesian" 

1060 FOR p=l TO 6 

1061 LET kp=K C p ) 

1062 IF p> =3 THEN LET kp=kp*RTD 

1063 PRINT kp ; TAB 16;XCp) 

1064 LPRINT kp 5 TAB l6;XCp) 

1066 NEXT p 

1068 PRINT “KEPLER-MEAN” ; TAB 1 6 ; "OSCULATING” 

1069 LPRINT "KEPLER-MEAN” ; TAB 1 6 ; " OSCULAT I NG ” 

1070 PRINT ADPJTAB 1 6 ; ao ’ EDP ; TAB 1 6 ; eo ’ I DP*RTD; TAB 1 6 ; i o*RTD ’ GDP 
*RTD; TAB 16;go*RTD’HDP*RTD;TAB 1 6 ; ho*RTD ’ LDP*RTD ; TAB 16;lo*RTD 

1071 LPRINT ADP ; TAB 1 6 ; ao ’ EDP ; TAB 1 6 ; eo ' I DP*RTD ; TAB 16;io*RTD'GD 
P*RTD; TAB 16;go*RTD’HDP*RTD»TAB 1 6 1 ho*RTD ’ LDP*RTD» TAB 16;lo*RTD 

1075 PRINT TAB 11; ”0— I ” : FOR p= 1 TO 6: PRINT TAB 7;DCp): NEXT p 

1076 LPRINT TAB 11;”0-I": FOR p= 1 TO 6: LPRINT TAB 7;DCp): NEXT 


p 

1080 LET DEL=SQR C DC 1 ) *DC 1 ) +D C 2 ) *DC 2 ) +DC 3 ) *DC 3 ) ) : IF DEL<=0.5E-3 
THEN GO TO 1095 

1085 IF DEL> =OLDEL THEN LET DEL=OLDEL : GO TO 1090 

1086 LET OADP= ADP : LET OEDP=EDP: LET OIDP=IDP: LET OGDP=GDP: LET 
OHDP=HDP : LET OLDP=LDP 

1087 LET OLDEL=DEL 

m 

1090 LET ADP=OADP : LET EDP=OEDP: LET IDP=OIDP: LET GDP=OGDP: LET 
HDP = OHDP : LET LDP = OLDP 
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1095 LET a0= ADP : LET eO=EDP: LET iO=IDP: LET g0=GDP: LET hO=HDP: 
LET 1 0=LDP 

1096 GO SUB 200 

1097 GO SUB 300 

1098 GO SUB 400 

1099 RETURN 

4000 REM [Converged Elements! 

4001 CLS : PRINT "[OSCULATING TO MEAN CONVERSION!" 

4002 LPRINT : LPRINT "[OSCULATING TO MEAN CONVERSION!" 

4005 PRINT " C Km-Km / sec-Deg . 3" 

4006 LPRINT " [ Km-Km / sec-Deg . 1" 

4010 PRINT "[MEAN ELEMENTS!";" DELR = " ; FN I(DEL);"km" 

4011 LPRINT "[MEAN ELEMENTS!";" DELR=";FN I (DEL); "km" 

4015 PRINT "a" = " ; ADP' "e ' * = M ;EDP' M i ' ' = " ; I DP*RTD * " g ' ' = " ; GDP*RTD ' " 
h * * = " ; HDP*RTD ' " 1 * ' ;LDP*RTD' ' 

4016 LPRINT "a' ' = " ; ADP ' "e ' ' = " ;EDP' “ i * ' =" ; I DP*RTD ' "g ' ' = “ ; GDP*RTD 7 
" h ' ' = " ; HDP*RTD ’ "1' ' = " ; LDP*RTD ' ' 

4017 PRINT "[OSCULATING ELEMENTS!" 

4018 LPRINT "[OSCULATING ELEMENTS!" 

4020 LET Ha=ADP*( 1 +e )-Re : LET Hp=ADP*( l-e)-Re 
4025 LET PD=2*PI /SQR ( GM / ADP A 3 ) : LET PD = PD/60 

4030 PRINT "x = ";x;TAB 1 6 ; " a= " ; a ' " y = " ; y ; TAB 1 6 ; " e = " ; e ' " z = " ; z ; TAB 
16; " i = " ; i *RTD ' "Dx = " ;Dx;TAB 16 ; M g=“ ; g*RTD * ”Dy = " ;Dy; TAB 16 J H h=" ;h* 
RTD' ,, Dz=“ ;Dz;TAB 16; "l = ";l*RTD 

4031 LPRINT " x=" ; x ; TAB 1 6 ; "a= " ; a ' " y = " ; y ; TAB 16; "e=" ;e ' "z=" ;z;TAB 
16; M i = " ; i*RTD' "Dx=" ;Dx;TAB 16; "g=" ;g*RTD' "Dy = " ; Dy ; TAB 16; M h=" ;h 

*RTD 7 " Dz = " ; Dz ; TAB 16; "1 = " ; 1 *RTD 

4035 PRINT : PRINT " r = " ; r ; TAB 1 6 ; " EA= " ; EA*RTD ' " V= " ; SQR V2;TAB 16 
; " f = " ; f *RTD 

4036 LPRINT : LPRINT "r= " ; r ; TAB 16; " EA= " ; EA*RTD * " V= " ; SQR V2;TAB 
16; "f = " ; f *RTD 

4040 PRINT "Hp=";Hp;TAB 16;" PD= " ; PD ; " m i n . " ' " Ha= “ ; Ha 

4041 LPRINT "Hp=" ;Hp; TAB 16;" PD= " ; PD ; " m i n . " ' " Ha= " ; Ha 
4999 RETURN 

9900 REM CONSTANTS 

9901 RESTORE 

9910 READ GM,Re,We , IDF, J2, J3, J4 

9920 DATA 398600.63,6378. 166,0.72921 1 59E- 4 , 298 . 25 , -0 . 10826517E-2 
,0. 25450306E-5,0. 16714987E-5 

9922 READ F 1 D2 , F 1 D3 , F 1 D4 , F3D2 , F3D8 , F3D32 , F5D4 , F5D 1 6 , F 1 5D 1 6 
9924 DATA . 5 , . 3333333333 , . 25 , 1 . 5 , . 375 , . 09375 , l . 25 , . 3 1 25 , . 9375 
9926 READ N 1 , N 2 , N3 , N 4 , N5 , N6 , N9 , N 1 1 , N25 , N30 , N35 , N40 , N90 , N 1 26 , N38 5 
9928 DATA 1,2,3,4,5,6,9,11,25,30,35,40,90,126,385 
9930 LET Ke =SQR <GM/Re^3) 

9940 LET K2=-0.5*J2*Re~2 
9950 LET K3 = J3*Re A 3 
9960 LET K4 = F3D8*J4*Re' v 4 

9980 LET RTD=180/PI: LET DTR=PI/180: LET f t Tkm=0 . 0003048 

9985 LET We =We *RTD 

9986 LET Rf=l/IDF: LET Zf=2*Rf- Rf*Rf 

9987 LET PI2=2*PI 

9989 RETURN 

9990 REM LOAD OSC ARRAY 

9992 DIM P< 2, 101,6) 

9993 DIM G ( 1 0 1 , 6 ) : DIM X(6): DIM K(6): DIM D<6): DIM 0(6) 

9994 RETURN 

9995 PRINT " BgOBSMEAN 

9996 PRINT "LOAD "ORBIT DATA" DATA PO" 

9997 LOAD "ORBIT DATA" DATA PC ) 

9998 LET DEL t = P ( 2 , 1 0 1 , 1 ) : LET END t = P ( 2 , 1 0 1 , 2 ) 

9999 RETURN 
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OSCULATING TO MEAN CONVERSION 

[ INPUT 1 


a = 6 7 4 3 . 68 2 


e = . 000107 


i = 28 .46 17 17 


g = 4 . 6 38 


h : - 2 5 5 . 318 


1 = 226 . 18 


I te rat i o 

n #: 1 

Kepler input 

Cartes i an 

6743.682 

-3365 . 0948 

.000107 

5287 .2238 

28.461717 

-2491 .0616 

4 .638 

-5.6410597 

255. 318 

-4.6817001 

226 . 18 

-2.31 48737 

KEPLER-MEAN 

OSCULATING 

6743.682 

6743.2323 

. 000107 

.00044416665 

28.461717 

28 . 458226 

4.638 

210.98579 

255 .318 

255 .35392 

226. 18 

19.8 14485 

o-r 


2 . 3029089 

5 .6256828 

-1 . 3741236 

.00047568977 

.0054707 1 17 

.001 1872994 

I terat ion # : 2 

KEPLER-MEAN 

OSCULATING 

6744 . 1 379 

6743.6818 

.00064203994 

.00010674693 

28. 465178 

28.461668 

22 . 57365 

5. 1222176 

255 . 28209 

255. 31806 

208 .26196 

225.69567 


0-1 

. 0076503754 
-.00051879883 
-.0034160614 
- 1 . 32 43 407E- 6 
-3. 3546239E-6 
-1 . 6335 398E-6 


I terat 
KEPLER-MEAN 
6744. 1381 
. 00064200723 
28 . 465227 
22 . 489951 
255 . 28203 
208 . 34577 

O-I 


ion # : 3 

OSCULATING 
6743.682 
.00010700 1 28 
28.461717 
4.6384579 
255 .318 
226 . 17954 


-.000026702881 

- .000053405762 

- .000014305 1 1 5 
9.31 32257E-9 
-6. 1 46729E-8 
2. 3283064E-8 


[OSCULATING TO MEAN CONVERSION! 

[ Kra-Km/sec-Deg . 1 

[MEAN ELEMENTS] DELR= . 00006 l 4km 
a’ ’=6744. 1381 
e ’ ' = . 00064200723 
i ' ' =28 . 465227 
g"=22. 489951 
h’ ’=255.28203 
1 ’ ' =208 .34577 


[OSCULATING ELEMENTS! 


x=-3365 . 0948 
y = 5 287 . 2238 
z=-249 1.0616 
Dx=-5. 6410597 
Dy=-4. 6817001 
Dz=-2 . 3148737 

r = 67 44 . 1817 
V=7. 6875556 
Hp=365. 25048 
Ha=366. 69373 


a = 67 43.682 
e= .00010699994 
i=28. 461717 
g=4. 6377326 
h = 255 . 3 1 8 
1=226. 18027 

EA=226 . 17584 
f = 226 . 17142 
PD=9 1 . 864838m i n . 
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OSCULATING TO MEAN CONVERSION 
t INPUT ] 
x=-0. 52760348 
y = 0. 82895001 
z = -0. 39056447 
Dx=-0. 71357062 
Dy = -0 . 59222684 
Dz = -0 . 2928 l 956 


I terat ion # : 1 
Kepler input 


6743.682 
. 00010673575 

28.461715 
4.6352318 

255. 31789 
226. 18342 
KEPLER-MEAN 
6743.682 
. 00010673575 

28.461715 
4.6352318 

255 .31789 
226 . 18342 


Cartes i an 
-3365. 1426 
5287 . 1808 
-249 1 . 085 
-5.6410252 
-4.6817602 
-2.314841 
OSCULATING 
6743 .2323 
.00044439565 
28 . 458224 
210.97213 
255 . 3538 1 
19.828806 


2.3028898 
5.6256752 
-1 . 3741741 
.00047560222 
.0054707192 
.001 1872184 


I tera 

KEPLER-MEAN 

6744. 1379 

.0006417772 

28.465177 

22.580998 

255.28198 

208 . 25527 


ion # : 2 

OSCULATING 

6743.6818 

.0001064839 

28. 461667 

5.1211681 

255.31795 

225.69738 


. 0075588226 
- .00056266785 
- .0034027 1 
-1.281 4999E-6 
-3. 4347 177E-6 
-1 . 5832484E-6 


I terat ion # : 3 


KEPLER-MEAN 
6744. 1382 
.00064174196 
28 . 465226 
22 . 497141 
255.28192 
208 . 33924 


OSCULATING 
6743.6821 
. 00010673641 

28. 461716 
4.6351024 

255. 31789 
226 . 18355 


O-I 

.000049591064 
-9 . 5367 432E-6 
3.81 46973E-6 
-3.7 25 2903E-8 
1 . 6763806E-8 
-1 . 676 3806E-8 


[OSCULATING TO MEAN CONVERSION! 

[ Km-Km / sec-Deg . 1 

[MEAN ELEMENTS! DELR= . 00005064km 
a' '=6744. 1382 
e' '=.00064174196 
i ' ' =28 . 465226 
g' '=22.497141 
h' '=255.28192 
1 ' ' =208 . 33924 


[OSCULATING ELEMENTS! 


x = - 3365 . 1426 
y = 5 287 . 1808 
z = - 2 49 1 .085 
Dx=-5. 6410252 
Dy=-4. 6817602 
Dz = - 2 . 314841 

r = 67 44 . 1805 
V=7. 687557 
Hp= 365 . 25237 
Ha=366 .69205 


a = 67 43.6821 
e=. 00010673623 
1=28.461715 
g=4. 63533 
h=255. 31789 
1=226. 18333 

EA = 2 26 . 1789 1 
f =226. 1745 
PD=9 1 .86484 
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7 

8 

100 
1 10 
1 15 
1 20 
1 25 

135 
1 40 
160 

199 

200 
201 

202 

203 

204 

205 

206 
207 


208 

209 


210 


OREM SAVE “ BqORBI T" 

2 REM An Economical Se m i - Ana 1 y t i ca 1 Orbit 

for Micro-Computer Applications 1 986-T 1 meX/S inc lair 0< 68 
R.A. Gordon, NASA Goddard Space Flight Center, GreenBelt, 

AIAA 24th Aerospace Sciences Meeting Jan . 6 9,1986 
[AIAA-86-00851 

4 DEF FN 1(S»C,T)=(PI AND C<0)+(PI2 AND C>0 AND S<0) 

+ CATN T AND COO) + (PI/2 AND C = 0 AND S>0) 

+ ( 3*PI / 2 AND C = 0 AND S<0) 

DEF FN M(X,N)=X-N*INT (X/N) 

DEF FN A( A ) =FN M(#i,PI2) n . 

DEF FN T(S,C)=(PI2 AND S<0)+(-l AND S<0 OR +1 AND S>-0) 

*ACS C 

REM SOLVE KEPLERS EQ . 

LET EA=0 

IF 1=0 THEN GO TO 160 

LET EA= 1 +e , 

FOR N=l TO 10: LET 0EA=EA : LET FE=EA-e*SIN EA-1: 

LET EA=EA-FE/( 1 -e *C0S ( EA-0 . 5*FE ) ) : LET DEA-ABS (EA-OEA) 

IF DEA< = 0 . 1 E- 8 THEN GO TO 160 
NEXT N 

LET EA=FN M(EA,2*PI> 

RETURN 

REM BgORBIT 
LET ADP=aO : LET EDP=eO: 

LET HDP=hO : LET LDP=10 
LET N0=SQR CGM/ADP^NS) 

LET EDP2=EDP*EDP : LET CN2=Nl-EDP2: LET CN=SQR 
LET GM2=K2 / ADP~N2 : LET GMP2=GM2 / ( CN2*CN2 ) : 

LET GM4=K4/ADP"N4: LET GMP4=GM4/CN~8 : LET F l D4G2=F 1 D4*GMP2 
IF Dt=0 THEN LET CI=C0S (IDP): LET CI2=CI*CI: 

LET Cl 3=CI 2*CI : LET CI4=CI2*CI2 
REM 1 DOT , gDOT , hDOT 

LET 1DOT=NO*(CN*(GMP2* (F3D2*(N3*CI2-N1) +F3D32*GMP2 

*(N25*CN2+16*CN-15+(N30-96*CN-N90*CN2)*CI2 

+(N25*CN2+144*CN+105)*CI4) ) )+F15D 1 6*GMP4*CN 
*EDP2*(N3-N30*CI 2+N35*CI 4) ) 

LET nDOT = NO + 1 DOT : LET PD = PI2/nD0T 
LET gDOT = NO* C ( GMP2*( F3D2*( N5*CI 2-N1 )+F3D32*GMP2 

*(N25*CN2+24*CN-N35+(N90-192*CN-126*CN2)*CI2 

+(45*CN2+360*CN+N385)*CI4) ) >+F5D16*GMP4 
*( 2 1-N9*CN2 + (N 1 26*CN2-270)*CI2 + (N385- 189*CN2)*CI 4) ) 

LET hDOT=NO*( ( GMP2* ( F3D8*GMP2* ( (N9*CN2+ 1 2*CN-N5> 
*CI-(N5*CN2+36*CN+N35)*CI 3)-N3*CI ) ) +F5D4*GMP4* ( N5-N3*CN2 ) 
*CI*(N3-7*CI2) ) 


LET I DP= i 0 : LET GDP=gO: 


( CN2 ) 


211 RETURN 

220 REM SECULAR-GDP , HDP , LDP 

221 LET GDP=gO+gDOT*Dt 

222 LET GDP=FN A (GDP ) 

223 LET HDP=hO+hDOT*Dt 

224 LET HDP=FN A ( HDP ) 

225 LET LDP=10+nD0T*Dt 

226 LET LDP=FN A(LDP) 

227 LET a = ADP : LET e =EDP : LET i = I DP : 

LET g=GDP : LET h=HDP : LET 1=LDP 

229 RETURN 

230 REM SP , LP- CONSTANTS 

232 LET CN3=CN2*CN : LET CN6=CN3*CN3: LET F 1 D 1 CN- 1 / < 1 +CN ) : 

LET F 1 QCN3= 1 /CN3 : LET F1DCN6=1/CN6 

233 LET GM3=K3/ADP~3: LET GMP3=GM3/CN6 : LET G3DG2=GMP3 /GMP2 

234 IF Dt=0 THEN LET SI=SIN (IDP): LET TI=SI/CI: 

LET P3T2M 1 -N3*CI 2-N 1 : LET P 1 MT2 = N 1 -Cl 2 : 

LET SQ 1MT2=SQR (P1MT2): LET T3 1 MT2 =N3*P 1 MT2 : 



LET T5T2M1=N5*CI2-N1 : LET P3M5T2=N3-N5*CI 2 : 

LET A0=CI2/(N1-N5*CI2) : 

LET A1=F1D2*F1D4*(N1-N1 1*CI2-N40*CI 2*A0) : 

LET A3=-F1D2*F1D4*CI*(N1 1 +80*AO+200*AQ*AO) 

235 LET EDPT3=N3*EDP : LET SP3=F 1 D2*GMP2 : LET TSP3=CI*SP3: 

LET SP6 = C I *SP3*SQ 1 MT2 

236 LET A2=CN3*GMP2*A1-F1D4*F1D4G2 

*( N2+EDP2- 400*EDP2*CI 2* AO* AO - 40* C N5*EDP2+N2 ) 

*CI 2*A0- 1 1*CI2*(N3*EDP2+N2) ) : 

LET A4=F t D4*G3DG2*S I : LET A5=( A4*EDP*CI ) / (N 1 +CI ) 

239 RETURN 

240 REM UDP, PERIODIC TERMS 

241 LET EP=EDP : LET GP=GDP: LET LP=LDP: LET UDP=GDP+LDP: 

LET UDP=FN A ( UDP ) 

242 REM LP-TERMS 

243 LET SG=SI N CGDP): LET CG=COS (GDP): 

LET S2G=N2*SG*CG : LET C2G=N2*CG*CG-N 1 

244 LET DIE=A4*SG+EDP*GMP2*A1*C2G: LET D 1 I =- < EDP*D 1 E ) /TI : 

LET D 1 E=CN2*D 1 E : LET D2E=EDP*CN3*GMP2*A 1 *S2G-CN3*A4*CG 

245 LET EP=SQR (D2E*D2E+( EDP + D 1 E) * < EDP+D 1 E) > 

246 LET HP=HDP+EDP2*A3*GMP2*S2G+( ( EDP*CI *A4)/(SI*SI ) ) *CG : 

LET HP=FN A ( HP ) 

247 LET UP=UDP+A2*S2G+( ( EDP*A4*F 1 D 1 CN ) * ( N2+CN-EDP2 ) +A5 ) *CG : 

LET UP=FN ACUP) 

248 LET SL=SI N (LDP): LET CL=COS (LDP) 

249 IF EDP> =0 . 05 THEN LET SM=D2E*CL+ ( EDP+D 1 E ) *SL : 

LET CM= ( EDP+D IE) *CL- (D2E*SL) : 

IF CMOO THEN LET TM=SM/CM: LET LP = FN 1(SM,CM,TM>: 

LET GP=UP-LP : LET GP=FN A(GP) : 

LET SG=SI N (GP): LET CG=C0S (GP) : 

LET S2G=N2*SG*CG : LET C2G=CG*CG-Nl 

250 REM FP 

251 LET 1 =LP : LET e=EP: GO SUB 100: LET EAP=EA: 

LET SEA=SI N (EA): LET CEA=COS (EA) 

252 LET ADR=N 1 / (N 1 -EP*CEA ) : LET ADR2=ADR*ADR : 

LET ADR 3= ADR 2* ADR : LET SF=ADR*SQR ( N 1 -EP*EP ) *SEA : 

LET CF=ADR* ( CEA-EP ) : LET FP=FN T(SF,CF) 

253 REM SP-TERMS 

254 LET CF2=CF*CF : LET CF3=CF2*CF: LET S2F=N2*SF*CF : 

LET C2F=N2*CF2-N 1 : LET S3F=N 3*SF-N4*SF*SF*SF : 

LET C3F=N4*CF3-N3*CF : LET S2GPF=S2G*CF+C2G*SF : 

LET S2GP2F=S2G*C2F+C2G*S2F : LET S2GP3F=S2G*C3F+C2G*S3F : 
LET C2GPF=C2G*CF-S2G*SF: LET C2GP2F=C2G*C2F-S2G*S2F : 

LET C2GP3F=C2G*C3F-S2G*S3F 

255 REM COMPUTE a,e,i,g,h,l 

256 LET a = ADP*( N 1 +GM2*( P3T2M 1 *< ADR3-F 1 DCN3 ) 

+T3 1 MT2*ADR3*C2GP2F ) ) 

257 LET D 1 E= ( F 1 D2*CN2* ( ( N3*F 1 DCN6*GM2*P IMT2*C2GP2F 
* ( EDPT3*CF2+N3*CF+EDP2*CF3+EDP ) ) - ( GMP2*P 1 MT2 
*(N3*C2GPF+C2GP3F) )+P3T2M 1 *GM2*F 1 DCN6*( EDP*CN+EDP 
*F1D1CN+EDPT3*CF2+N3*CF+EDP2*CF3) ) )+DlE: 

LET D2E=-F1D4G2*CN3*(N2*P3T2M1 * ( ADR2*CN 2+ADR+N 1 )*SF+T31MT2 
*( (-ADR2*CN2-ADR+N1 ) *S2GPF+ ( ADR2*CN2+ADR+F 1D3 ) *S2GP3F ) ) 
+D2E : LET e =SQR ( D2E*D2E+ (EDP+D 1 E) *(EDP+D IE) ) 

258 LET i = I DP+D 1 I +SP6*(N3*C2GP2F + EDPT3*C2GPF+EDP*C2GP3F) : 

LET i =FN A ( i ) 

259 LET h=HP-TSP3* ( N6 * ( FP-LP + EDP*SF ) 

- ( N 3*S2GP2F+EDPT3*S2GPF+EDP*S2GP3F ) ) : LET h=FN A(h) 

260 LET u = UP‘+( F 1 D 1 CN*F 1 D4G2*EDP*CN2* ( T31MT2*(S2GP3F*(F1D3 

+ADR2*CN2+ADR ) +S2GPF* ( N l - ( ADR2*CN2+ADR ) ) ) +N2*SF 
*P3T2M1*( ADR 2 *CN 2+ADR+N 1 ) ) ) +GMP2*F3D2* ( T5T2M 1 
*(EDP*SF+FP-LP) )+P3M5T2*(FID4G2*(EDP*S2GP3F+N3 
*(S2GP2F+EDP*S2GPF) ) ) : LET u=FN A(u) 
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261 LET SM=D2E*CL+(EDPtDlE)*SL: LET CM=(EDP+DIE) *CL-D2E*SL : 

IF CM< > 0 THEN LET TM = SM/CM 

262 LET 1=FN l < SM , CM , TM ) 

264 LET g = u-l : LET g=FN A(g) 

269 RETURN 

270 REM ORBGEN 

27 i LET Dt = t : IF Dt = 0 THEN GO SUB 200: GO SUB 230 
2 7 2* IF DQTa = 0 THEN GO TO 290 

273 IF t=0 THEN LET tsum=0: LET t0=0: GO TO 289 

274 LET s i gn= 1 : LET Dt*t-t0: IF Dt<0 THEN LET sign=-l 

275 LET PD=sign*PD - 

276 LET nPD= 1 : IF ABS Dt > ~ ABS PD THEN LET nPD= I NT (Dt/PD) 

277 IF ABS t < ABS (tsum + PD) THEN GO TO 289 

278 FOR n ~ 1 TO nPD 

279 LET DQTe = ( < 1 -eO) /aO)*DOTa 

280 LET D0TnD2= -(3/4)*<nQ/aO) *DOTa 

281 LET DELa=DOTa*PD : LET DELe =DOTe *PD : LET DELI =D0TnD2*PD*PD 

282 LET aO=aO+DELa : LET eO=eO+DELe: LET 10=10+DEL1: 

LET 1 0 = FN M ( 1 0 , 2*P I ) 

283 LET Dt=PD: LET tsum=tsum+PD 

284 GO SUB 200: GO SUB 220 

285 LET a0=a : LET eO-e: LET gO = g: LET hO=h: LET 10=1 

287 NEXT n 

288 GO SUB 200: GO SUB 230 

289 LET 1 0 = t s u m : LET Dt = t-tO 

290 GO SUB 220: GO SUB 240 
299 RETURN 

9900 REM “ CONSTANTS" 

9901 RESTORE 

9910 READ GM,Re , We , IDF, J2 , J3, J4 

9920 DATA 398600 . 6 3 , 6 378 . 1 66 , 0 . 7 29 2 1 1 59E- 4 , 298 . 25 , 

-0 * 108265 17E- 2, 0. 25 450306E-5 , 0 . 167 14987E-5 
9922 READ F 1 D2 , F 1 D3 , F 1 D 4 , F3D2 , F3D8 , F3D32 , F5D4 , F5D 1 6 , F 1 5D 1 6 
9924 DATA . 5 3333333333, . 25 , 1 * 5 375 09375 , 1 . 25 3 1 25 9375 
9926 READ N1,N2,N3,N4,N5,N6,N9,N1 l ,N25,N30,N35,N40,N90,N126,N385 
9928 DATA 1,2,3,4,5,6,9,11,25,30,35,40,90,126,385 
9930 LET Ke =SQR ( GM / Re A 3 ) 

9940 LET K2=-0*5*J2*Re*2 
9950 LET K3 = J 3*Re A 3 
9960 LET K4=F3D8*J4*Re~4 
9970 LET PI 2= 2*P I 

9980 LET RTD= 1 8 0 / P I : LET DTR = PI/180 
9999 RETURN 


D-3 


State matr i zar t-OC J2 ) 


300>REM Secular retarded 

301 DIM SC7,7> 

302 LET p=a0*cn2; LET p2=p*p: LET B2=C2/p2 

303 LET lmDOT=B2*cn*<3*CI2-l > *n0 

304 LET gmD0T = B2*C 5*CI 2 - 1 ) *n0 

305 LET hmDOT=- 2*B2*CI *n0 

306 LET FlDA=l/a0: LET FEDP=eO/p 

307 LET S 1 =F7D2*F 1 DA : LET S2=a0*FEDP: LET S3=2*B2*nO*SIN iO: 

LET S4=F3D2*FlDA*n0 

308 LET S( 1 , 1 ) = 1 

309 LET SC 2 , 2 ) = 1 

310 LET S C 3 , 3 ) = 1 

31 1 LET SC 4, 4) = i 

312 LET S C 5 , 5 ) = 1 

313 LET SC 6,6) = 1 

319 RETURN 

320 REM t 

321 REM g partials 

322 LET S(4,l)= -Sl*gmDOT*t 

323 LET S ( 4 , 2 ) = 4*S2*gmDOT*t 

324 LET S C 4 , 3 ) = -5*S3*CI*t 

325 REM h partiaLs 

326 LET SC 5 , 1 ) =- Sl*hmDOT*t 

327 LET S C 5 , 2 ) = 4*S2*hmD0T*t 

328 LET SC 5 , 3) = S3*t 

329 REM 1 partials 

330 LET SC6, 1 ) =- C S4+S 1 *1 mDOT ) *t 

331 LET SC 6 , 2 ) = 3*S2*lmDOT*t 

332 LET SC6>3)= - 3*cn*S3*CI * t 
340 RETURN 

350 REM DOTa 

352 LET F 1 DA= 1 /aO : LET pO=aO*C 1 -e0*e0) : 

LET FEDP=eO/pO: LET DNDA=nO/aO 

353 LET D4=Cl-eO): LET D3=FEDP*D4: LET D2=D0Ta*e0: 

LET D 1 =D0Ta*D3*F 1 DA 
355 LET S C 7 , 7 ) = 1 

369 RETURN 

370 REM retarded partials 

371 LET 1 2 = t*t 

372 LET SC 1,7)= t 

373 LET SC 2 , 1 ) = -D4*F 1 DA*F 1 DA*DOTa* t 

374 LET S C 2 , 2 ) = 1 -DOTa*F 1 DA*t 

375 LET S C 2 , 7 ) = D4*FlDA*t 

378 LET SC 4, 1 )=SC4, 1 ) - 4*D 1 *gmDOT*t 2 

379 LET SC4,2)=SC4,2)+2*D2*gmD0T*t2 

380 LET SC4,7)=SC4, 1 ) *t+4*D3*gmD0T*t 2 

382 LET SC5, i)=SC5, 1 ) -4*D1 *hmDOT*t 2 

383 LET SC5, 2 ) =S C 5 , 2 ) + 2*D2*hmDOT* t 2 

384 LET SC5,7)=SC5, 1 ) * t +4*D3*hmD0T*t 2 

386 LET SC6, 1)=SC6, 1)+F5D4*F3D2*F1 DA*DNDA*DOTa*t 2 

387 LET SC6, 2)=SC6, 2)-D2*lmDOT 

388 LET SC6,7)=SC6, 1 ) *t-F3D2*F 1 D2*DNDA*t2+3*D3*l mDOT*t 2 
399 RETURN 
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